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摘 要 


本 文 研究 了 分 数 微分 方程 和 非 线性 发 展 方程 中 的 若干 求解 的 方法 ， 

1. 高 维 Mittag-Leffler 函数 的 构造 及 其 积分 福 质 的 研究 ; 

2. 分 数 微 分 形式 空间 的 研究 及 分 数 梯 度 ， 旋 度 和 和 散 度 的 构造 ， 

3. 一 种 Darboux 变换 的 构造 及 其 在 非 线 性 发 展 方程 中 的 应 用 . 

第 一 章 介 绍 了 论文 的 基础 知识 以 及 前 人 工作 的 成 果 : 包括 介绍 数学 机 械 化 的 思想 与 
应 用 的 情况 ， 回 顾 分 数 微 积分 的 历史 和 所 应 用 的 领域 ， 简 单 介绍 孤立 子 研究 的 历史 、 发 
展 与 成 果 ， 

第 二 章 主 要 研究 了 分 数 微分 方程 和 分 数 微分 形式 的 一 些 问 题 ， 

(1) 介绍 了 分 数 微分 方程 研究 的 历史 发 展 和 现状 并 且 构 造 了 高 维 空间 中 的 Mittag- 
Lefer BR, 研究 了 它 的 一 些 重要 的 积分 变换 , 包括 Laplace 变换 ，Mellin WR, W Fourier 
变换 ， 并 且 用 它 来 求 分 数 扩散 - 波动 方程 的 基本 解 ， 还 与 经 典 的 扩散 和 波动 方程 进行 了 
XE. 

(2) 介绍 了 分 数 微分 形式 的 理论 知识 和 发 展 , 并 将 分 数 微分 形式 和 经 典 的 微分 形式 作 
比较 . 进一步 ， 为 了 物理 学 上 的 需要 ， 试 图 用 分 数 外 导数 在 分 数 微分 形式 空间 中 构造 出 
RRR, PERRE. 

第 三 章 给 出 了 “ AC-BD " 理论 的 基本 思想 和 应 用 , 以 及 构造 CD 对 的 理论 与 算法 ， 
在 这 个 思想 的 指导 下 ， 我 们 构造 了 一 种 Darboux 变换 ， 并 利用 这 个 变换 给 出 了 (1+1)- 维 
高 阶 Broer-Kaup(HBK) 系统 的 一 些 新 解 . i 


关键 词 : 数学 机 械 化 ; 分 数 微分 方程 ; 分 数 微分 形式 ，Mittag-Leffer BAK, Darboux 变换 
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Some Problems in Fractional Differential Equations and 
Nonlinear Evolution Equations 


Abstract 

This dissertation primarily considers the following problems of fractional differential equa- 
tions and nonlinear partial differential equations. 

1. Construction of Mittag-Leffler type function and studies on its integral properties. 

2. Some studies on the Fractional differential forms and of fractional gradient, curl and 
divergence. 

3. Introduction of a kind of Darboux transform, and verification of its correctness by virtue 
of the Computational Symbolic Software Maple and display of its applications in the nonlinear 
evolution equations. 

Chapter one involves a survey of the fundamental knowledge and achievement obtained by 
foregoers on the subjects for the purpose of the later parts in this dissertation, including the 
introduction of the history and development of mathematics mechanization, brief review of the 
history and status in quo of the fractional calculus and soliton theory. 

Chapter two focuses on several problems of fractional differential equations and fractional 
differential forms: 

(1) Constructs a new type of Mittag-Leffler function, gives some of its important integral 
transforms such as Laplace transform, Mellin transform and Inverse-Fourier transform. And 
then emphasizes on its application on the fractional diffusion-wave equations. 

(2) Based on the ideas and work of fractional exterior derivative and fractional differential 
forms, compares the fractional differential forms them the standard ones. Furthermore, seeks 
the fractional gradient, curl and divergence in the fractional differential form space in order to 
meet the develop of theory in physics and finance. 

Chapter three of this dissertation is devoted to introducing the essentials of "AC=BD”. 
Also, under the guidance of such idea we try to construct a kind of Darboux transform and 


apply to the (1--1)-dimensional higher Broer-Kaup(HBK) system. 


Keywords: Mathematics Mechanization; Fractional differential equations; Fractional differen- 
tial forms; Mittag-Leffler function; Darboux function. 
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独创 性 说 明 


作者 郑重 声明 : 本 硕士 学 位 论文 是 我 个 人 在 导师 指导 下 进行 的 
研究 工作 及 取得 研究 成 果 . 尽 我 所 知 ,除了 文中 特别 加 以 标注 和 至 
澳 的 地 方 外 ,论文 中 不 包含 其 他 人 已 经 发 表 或 撰写 的 研究 成 果 ， 也 
不 包含 为 获得 大 连理 工大 学 或 其 他 单位 的 学 位 或 证 书 所 使 用 过 的 材 
料 。 与 我 一 同 工 作 的 同志 对 本 研究 所 做 的 贡献 均 已 在 论文 中 做 了 明 
确 的 说 明 并 表示 了 谢意 . 
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大 连理 工大 学 学 位 论文 版 权 使 用 授权 书 


本 学 位 论文 作者 及 指导 教师 完全 了 解 “ 大 连理 工大 学 硕士 、 博 士 学 位 论文 版 权 使 用 
规定 "， 同 意 大 连 理工 大 学 保留 并 向 国家 有 关 部 门 或 机 构 送 交 学 位 论文 的 复印 件 和 电子 
版 ， 允 许 论 文 被 查阅 和 人 惜 阅 。 本 人 授权 大 连理 工大 学 可 以 将 本 学 位 论文 的 全 部 或 部 分 内 


容 编 入 有 关 数 据 库 进行 检索 ， 也 可 采用 影印 、 缩 纯 或 扫描 革 复 制 手段 保存 和 汇编 学 位 论 
文 。 


(2% 


作者 签名 : 


7 
导师 签名 : Ic 


40 


大 连理 工大 学 硕士 学 位 论文 


第 一 章 H i 
81.1 。 数学 机 械 化 与 符号 计算 


回顾 数学 的 发 展 历史 ， 主 要 有 两 种 思想 ， 一 个 是 公理 化 的 思想 ， 另 一 个 是 机 械 化 的 
思想 ， 前 者 源 于 希腊 ， 后 者 则 贯穿 整个 中 国 古 代数 学 ,这 两 种 思想 彼此 辉映 ,贯穿 整 个 
数学 历史 ， 对 数学 的 发 展 都 曾 起 到 了 巨大 的 作用 ， 从 汉 初 完成 的 《 九 章 算术 》 中 对 开平 
方 、 开 立方 的 机 械 化 过 程 的 描述 , 到 宋 、 元 时 代 发 展 起 来 的 求解 高 次 代数 方程 组 的 机 械 化 
算法 ， 无 一 不 与 数学 机 械 化 思想 有 关 ， 对 数学 的 发 展 起 到 了 巨大 的 作用 . 公理 化 的 思想 
在 现代 数学 ， 特 别 是 在 纯粹 数学 中 占据 着 统治 地 位 . 但 即使 在 现代 纯粹 数学 的 研究 中 ， 
机 械 化 的 思想 也 一 直 发 挥 着 重要 的 作用 ，Hilbert 所 倡导 的 数理 逻辑 为 计算 机 的 设计 原理 
做 了 准备 . 数学 巨匠 E. Cartan 关于 微分 方程 、 微 分 几何 及 李 群 的 著作 中 也 显现 了 机 械 化 
的 特色 。 H. Cartan 关于 代数 拓扑 学 中 同调 群 计算 的 工作 也 可 以 看 作 是 机 械 化 思想 成 功 
的 范例 ， 

在 以 计算 机 为 标志 的 信息 革命 时 代 到 来 之 际 ， 数 学 应 该 有 什么 样 的 创新 与 之 相 适 应 
E? 正 是 基于 这 种 考虑 , 中 国 科学 院 院士 吴 文俊 先生 倡导 数学 机 械 化 的 研究 | 21, 20 til 
纪 70 年 代 ， 吴 文俊 由 中 国 的 传统 思想 出 发 ， 从 几何 定理 证 明 入 手 开始 数学 机 械 化 研究 ， 
所 建立 的 数学 机 械 化 方法 ， 不 仅 将 中 国 传统 数学 发 扬 光 大 ， 而 且 也 为 国际 自动 推理 的 研 
究 开 辟 了 新 的 前 景 . 经 过 近 20 年 的 努 方 ， 几何 定理 自动 证 明 的 吴 方 法 及 在 其 影响 下 产生 
的 一 系列 重要 的 新 方法 ， 已 经 发 展 成 就 有 我 国 特色 在 国际 上 领先 的 数学 机 械 化 理论 ,这 
一 理论 已 不 仅 在 几何 定理 的 机 器 证 明 ， 方 程 组 求解 ， 微 分 几何 ， 理 论 物理 ， 力 学 等 领域 
得 到 成 功 应 用 ， 还 为 机 器 人 学 ,数控 技术 ， 几 何 辅助 设计 ，CAD, 计算 机 视觉 等 高 科技 领 
域 提供 了 有 力 工具 ,他 引入 的 求解 非 线性 代数 方程 组 的 吴 方 法 是 求解 代数 方程 组 精确 解 
”最 完整 的 方法 之 一 ， 已 经 被 成 功 地 用 于 解决 很 多 问题 ， 并 实现 在 当前 流行 的 符号 计算 软 
件 中 . 上 世纪 80 年 代 末 ， 在 RittB, 和 等 人 工作 的 基础 上 ， 吴 先生 进一步 给 出 了 吴 微 分 消 
元 法 ， 提 出 了 吴 微 分 特征 列 的 概念 ， 完 善 和 发 展 了 特征 集 理论 ， 

高 小 山 研究 员 、 张 景 中 院士 、 周 成 青 教授 9. 6 合作 提出 了 基于 几何 不 变量 的 “ 消 点 
法 *， 实 现 了 自动 生成 几何 定理 可 读 证 明 这 一 目标 . 依据 这 一 方法 编制 的 程序 已 经 用 于 证 
明了 数 百 条 定理 。 其 中 大 部 分 的 证 明 是 简短 可 读 的 。 这 一 方法 不 仅 提高 了 定理 证 明 的 质 
量 ， 从 而 在 理论 上 极 大 地 推进 了 自动 推理 的 研究 ， 还 被 用 来 解决 了 CAD, 智能 CAI( 计 算 
机 辅助 教学 ) 与 机 器 人 中 的 若干 关键 理论 问题 . 这 一 工作 被 认为 是 五 十 年 代 Gerlentner 的 
经 典 工作 及 七 十 年 代 吴 方法 以 来 这 一 领域 又 一 重要 进展 . 美国 数学 会 “自动 定理 证 明成 
就 奖 * R “J. Macarthey 程序 验证 奖 ” 得 主 Boyer 称 , 该 工作 “是 自 ( 五 十 年 代 ) Slagle 与 
Moses 符号 积分 程序 以 来 自动 推理 界 最 重要 的 一 件 单独 的 事情 ?. 该 工作 “在 使 计算 机 象 
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具有 算术 天 才 那 样 具 有 几何 天 才 这 一 不 可 避免 的 过 程 中 将 是 一 座 里 程 碑 *. 自动 推理 届 权 
威 Loveland 在 AI Magazine 的 文章 中 将 这 一 工作 列 为 近年 来 自动 推理 界 “ 重 要 进展 ” 的 
第 一 项 ， 称 “在 几何 中 证 明 有 意义 的 定理 并 给 出 可 读 证 明 ( Gelentner 五 十 年 代 的 重要 工 
作 ) 近 年 来 才 由 周 威 青 , 高 小 山 , 张 景 中 的 几何 定理 证 明 器 超过”. 在 几何 自动 推理 方面 ， 
他 们 提出 了 微分 几何 的 自动 定理 证 明 的 新 方法 ， 所 编程 序 自动 证 明了 上 百 个 定理 并 发 现 
THAR. 给 出 了 Caley-Klein 几何 的 转换 定理 大 大 简化 了 非 欧 几何 的 自动 定理 证 明 , 用 
自己 研制 的 程序 解决 了 Zassenhaus 与 Maclane83 公开 问题 ， 这 充分 说 明了 他 们 的 方法 的 
有 效 性 .提出 了 几何 推理 的 演绎 数据 库 方法 ， 改 进 了 基于 搜索 的 定理 证 明 方 法 ， 并 第 一 
次 用 此 类 方法 证 明了 大 量 几何 定理 , 吴 尽 昭 研究 员 、 刘 卓 军 研究 员 U^ 9 将 吴 代数 消 元 法 
运用 到 逻辑 中 去 ， 较 好 地 解决 了 逻辑 中 一 阶 定理 证 明 的 问题 . 石 替 研 究 员 10.11 利用 吴 
代数 方法 ,研究 了 著名 的 Yang-Baxter 方程 的 解 的 问题 . 之 后 ， 他 利用 张 鸿 庆 教授 提出 的 
* AC-BD "的 思想 ， 将 Yang-Mills 方程 约 化 为 三 个 简单 的 二 阶 线性 微分 方程 。 王 世 坤 
研究 员 、 吴 可 研究 员 (12, 15) 将 吴 方 法 应 用 于 研究 Yang-Baxter Hl ( 包括 带 参数 、 带 色 参 
数 、 带 谱 参数 等 ) 的 解 的 结构 问题 . 

ERRAR Y 1141-06] 利用 吴 代 数 消 元 法 ， 在 求解 孤立 子 方程 方面 做 了 很 多 的 出 色 
IF. 李 志 滤 教授 和 石 替 研 究 员 利用 吴 方 法 成 功 地 得 到 了 BZ 方程 的 6 种 行 波 解 ， 朱 思 
网 教授 071 根据 AMS 猜测 ， 利 用 符号 计算 和 吴 代 数 消 元 法 ， 对 偏 微分 方程 的 Painleve 性 
质 进行 了 研究 , 证 明了 一 批 方程 具有 P- ER. 范 轧 贵 教授 在 这 方面 也 做 了 大 量 的 工作 ， 
他 推广 了 tanh 方法 ,借助 于 计算 机 并 利用 吴 方 法 得 到 了 很 多 方程 的 精确 解 。 朝 各 教授 运 
用 吴 微分 特征 集 求解 Lie 对 称 并 获得 了 一 批 对 称 ， 并 给 出 了 算法 和 程序 ， 成 功 地 将 吴征 
分 的 方法 运用 在 偏 微分 方程 {组 ) 的 求解 中 . 

近年 来 ， 张 鸿 庆 教授 及 其 课题 组 成 员 在 微分 代数 方程 的 代数 化 和 机 械 化 方面 做 了 大 
量 的 工作 08 | 张 鸿 庆 教授 早 在 1978 年 ， 将 代数 化 的 思想 引入 了 微分 方程 中 ， 给 出 了 具 
体 的 代数 化 算法 ， 解 决 了 大 量力 学 中 方程 组 的 解 的 恰当 性 问题 。 随 后 ， 证 明了 非 齐 次 线 
性 算 子 方程 组 Au=f 的 一 般 解 为 u=Cv+e ， 其 中 v 满足 方程 组 Dv-g, D 是 对 角 和 矩阵 ， 
用 代数 方法 给 出 了 C., D, e 的 具体 构造 方法 . 用 此 方法 可 以 给 出 各 种 弹性 力学 位 移 函 
数 和 应 力 函数 的 机 械 化 算法 ， 构 造 数学 物理 中 一 系列 方程 的 一 般 解 ， 基 于 吴 微 分 消 元 理 
$, 给 出 了 “ AC=BD ”除法 ， 利 用 这 一 除法 ， 可 以 得 到 一 些 方程 的 变换 ， 进 而 通过 变换 
RHE. SW 419 ， 陈 勇 博 士 20 在 微分 方程 的 求解 代数 化 方面 做 了 大 量 的 工 
作 . 他 基于 两 种 Riccati 方程 ， 提 出 了 新 的 求解 非 线性 微分 方程 的 更 有 效 的 方法 ， 获 得 了 
丰富 的 精确 解 . 最 近 ， 我 们 课题 组 成 员 又 将 吴 方 法 与 “ AC-BD " 理论 相 结 合 ,编制 了 广 
X tanh 方法 求解 ， 广 义 Riccati 方程 展开 法 求 类 孤立 子 解 ， Lie 对 称 及 其 Painleve 检验 软 
件 包 . 郑 学 东 硕士 CH ， 完 成 了 三 角 函 数 方法 求 精确 的 软件 包 ， 广 义 tanh 函数 求解 的 软 
件 包 ， 和 Painleve 检验 的 程序 包 ， 吕 卓 生 博士 完成 了 Lie 对 称 的 程序 编制 . 
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计算 机 代数 (Computer algebra 或 Symbolic and algebra computation) 兴起 于 上 世纪 60 
年 代 初 , 是 介 于 计算 机 , 数学 与 人 工 智能 之 间 的 一 门 边缘 学 科 , 计算 机 代数 的 主要 研究 内 
容 是 计算 机 上 数学 公式 演绎 的 算法 和 系统 应 用 ， 其 发 展 ， 大 致 可 分 为 以 下 儿 个 阶段 221: 

1, 20 世纪 60 E: 发 展 初期 着重 于 多 项 式 算术 ， 积 分 方法 ; 

2. 20 世纪 70 年 代 ; 诞生 了 Macsyma, Reduce, 以 及 用 于 抽象 域 的 Scratchpad /I 

3. 20 世纪 80 年 代 : 多 项 式 时 间 方 法 ， 因 式 分 解 等 的 研究 ， Maple 与 Mathematica 相 
RRA, 

4. 20 世纪 90 年 代 ， 微 积分 ， 数 学 的 Web 发 布 ， 黑 箱 符号 对 象 

5. 本 世纪 初 :符号 计算 ， 数值 计算 ， 几何 ， 组 合 与 逻辑 范例 的 融合 . 

计算 机 代数 的 主要 功能 在 于 ， 在 计算 机 上 以 符号 形式 进行 运算 ， 实 现 公式 的 机 器 推 
R. 如 果 和 初等 数学 相 比 ， 寻 常 的 数值 计算 可 比拟 为 数值 计算 ， 而 符号 计算 可 比拟 为 代数 
运算 ( 故 有 计算 机 代数 之 说 ), 前 者 只 适合 于 个 例 分 析 , 后 者 的 结果 则 具有 普遍 性 . 由 于 计 
算 机 代替 了 人 工 的 推导 , 演算 速度 成 千 万 们 地 增加 , 使 得 原来 令 人 望 而 生 长 的 复杂 计算 变 
得 简易 快 洁 , 在 最 初 的 十 几 年 间 , 计算 机 代数 也 被 叫做 : 符号 与 代数 计算 、 公式 处 理 、 机 
器 代数 等 ,于 是 在 1979 年 ， 美 国 加 州 理 工学 院 的 教授 统一 了 计算 机 代数 的 名 称 和 定义 : 

“计算 机 代数 是 一 门 利用 数学 、 计 算 机 进行 代数 和 解析 处 理 或 操作 的 学 科 ”。 1982 年 ， 
德国 卡尔 斯 鲁 厄 大 学 的 RR. Loose 教授 又 进行 了 进一步 的 拓 广 和 补充 ，“ 计 算 机 代数 是 设 
it. 分 析 、 改 造 和 应 用 代数 算法 的 计算 机 科学 分 支 *， 计算机 代数 的 最 早出 现 公认 以 1960 
年 美国 麻 省 理工 学 院 的 Mccarthy 推出 的 Lisp 语言 为 标志 . 在 随后 的 十 几 年 间 ， 计 算 机 代 
数 的 发 展 引起 了 国际 计算 机 科学 界 的 重视 .美国 的 计算 机 协会 组 织 了 符号 与 代数 处 理 专 
业 组 ( Special interest Group on Symbolic and Algebraic Manipulation ， 简 称 SIGSAM ), 其 
RAGA 30 多 个 国家 .这 个 国际 组 织 每 两 年 召开 一 次 国际 会 议 ( 简称 ISSAC ), 专 
门 交流 计算 机 代数 方面 的 研究 成 果 ， 西欧 各 国 计 算 机 工作 者 组 织 了 欧洲 符号 和 代数 处 理 
专业 委员 会 ( Symbolic and algebraic Manipulation of European ， 简 称 SAME )， 定 期 召开 
简称 AAECO 的 国际 会 议 . 这 两 大 组 织 还 分 别 创 办 了 计算 机 代数 刊物 SIGAM Bulletin 和 
Journal of Symbolic computation .这 些 学 术 活 动 大 大 地 推动 了 计算 机 代数 的 发 展 . 

在 过 去 的 三 十 儿 年 间 , 计算 机 代数 取得 了 诸多 的 成 就 23] Risch PAIRT: 数学 函 
数 在 闭 形式 下 的 积分 问题 是 可 判定 的 。 Berlekamp 29 给 出 了 有 效 分 解 多 项 式 ( 模 一 个 大 
XO 的 随机 化 算法 . 因为 多 项 式 的 系数 域 的 不 同 ， 这 时 需要 给 出 一 个 更 一 般 的 算法 . 在 
上 个 世纪 的 七 十 年 代 的 早期 ， Berlekamp 给 出 了 这 个 代数 算法 在 抽象 代数 域 上 的 更 一 般 
的 程序 ，Brownl26] 给 出 了 新 的 有 效 的 多 元 多 项 式 最 大 公 因子 的 符号 计算 方法 . Gosper??7] 
对 于 无 限 超 几何 和 给 出 了 一 个 精巧 的 算法 。 Lenstra et al 29: 291 的 Lovase's 的 格 点 约 化 
算法 ， 是 具有 深远 意义 的 Euclidean 算法 的 一 般 化 ， 最 初 是 作为 一 个 多 项 式 分 解 的 子 算法 
出 现 的 , 稀 焉 多 元 多 项 式 的 差 补 算法 名 (其 中 一 部 分 是 基于 修正 误差 码 的 ) 已 经 变 成 了 在 
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计算 黑匣子 多 项 式 中 的 基本 工具 (Kaltofen & Trager), 1995 年 ， Shoupl51 在 修改 了 由 
Kaltofen & Shoup 给 出 的 关于 一 元 多 项 式 (系数 在 有 限 域 上 ) 的 分 解 算法 后 ， 使 这 个 算法 
更 加 实用 化 . 特别 是 在 上 个 世纪 的 90 年 代 ， 借 助 于 一 些 符号 计算 软件 工具 ， 计 算 机 代数 
的 发 展 更 加 迅猛 ， 而 且 不 断 地 应 用 到 其 他 领域 ， 如 高 能 物理 、 天 体力 学 、 广 义 相 对 论 、 电 
子 光学 、 分 子 物理 、 自 动 化 、 航 空 学 、 生 物 学 和 化 学 等 等 ， 


81.2 ”分 数 微 积 分 历史 和 发 展 概况 


数学 是 一 门 美丽 的 自然 语言 ， 但 它 却 常常 会 向 我 们 传递 一 些 误导 人 的 信息 ， 比 如 说 
我 们 常用 的 自然 数 和 实数 一 自然 数 其 实 并 不 自然 ， 绝 大 多 数 的 实数 在 现实 生活 的 计算 
中 也 并 不 是 真实 存在 的 分数 微 积分 这 个 概念 也 是 这 样 一 个 “名 不 符 实 ” 的 数学 名 词 : 
它 不 是 指 分 数 的 微 积分 ， 也 不 是 光 指 分 数 阶 的 微 积分 学 ， 而 是 代表 任意 阶 积分 和 微分 的 
名 词 ， 并 且 推 广 、 统一 了 整数 阶 

分 数 微 积分 的 发 展 历史 几乎 与 整数 阶 微 积分 的 发 展 史 相同 ， 自 从 分 数 导数 的 想法 最 
HL L'Hospital F 1695 年 给 G.W.Leibnitz 的 一 封 信 中 提出 ， 人 们 对 这 个 学 科 的 兴趣 就 从 
未 停止 过 .尽管 Euler Lagrange 和 其 它 的 数学 家 曾经 在 更 早 的 时 期 在 这 方面 做 出 了 贡 
献 ， 但 是 这 方面 最 早 的 系统 研究 则 是 由 Liouville Riemann 以 及 Holmgren 在 19 世纪 初 
叶 与 中 叶 完 成 的 。 最 早 是 Liouville 将 水 数 展开 成 指数 形式 ， 并 通过 逐 项 运算 定义 了 这 种 
展开 形式 的 q 阶 导数 ( 这 里 gq 是 正 整数 )。 Riemann 提出 了 一 种 不 同 的 定义 ( 对 给 定 的 
WH), 这 个 定义 可 以 应 用 于 带 有 负 短 次 项 的 震级 数 形式 , 最 后 是 Grunwald 和 Krug 第 一 
次 将 Liouville 和 Riemann 的 定义 统一 了 起 来 . 参见 92.33.34]. 

与 这 些 理论 的 起 始 阶段 平行 ， 对 各 种 各 样 的 问题 分 数 微 积分 的 应 用 取得 了 发 展 , 在 
一 定 意义 上 ， 分 数 微 积分 最 初 是 由 Abel 于 1823 年 发 现 的 ， 他 发 现 对 于 tautochrone ， 其 
积分 方程 的 解 可 以 通过 一 个 积分 变换 得 到 ， 正 如 我 们 将 要 看 到 的 那样 ， 这 主要 得 益 于 问 
十 被 写成 半 导 数 的 形式 ， 用 分 数 微 积分 解决 问题 的 一 个 强 有 力 的 触动 来 自 于 Boole 的 符 
号 方法 的 发 展 ， Boole 用 这 种 方法 来 求解 常 系数 的 线性 微分 方程 ， Boole 思想 的 本 质 是 
将 微分 算 子 展开 成 任意 函数 的 形式 第 级 数 ， 于 是 微分 方程 的 解 就 是 这 些 级 数 的 形式 道 . 
这 种 方法 对 某 些 特定 的 函数 类 的 要 求 特 别 严格 ， 而 且 已 经 在 许多 方面 都 得 到 了 推广 . 

然而 ， 由 于 长 期 没有 得 到 实际 应 用 背景 的 促进 而 发 展 缓慢 ， 直 到 Heaviside 发 展 了 算 
子 微 积分 ， 用 于 求解 某 些 电磁 理论 方面 的 问题 ， 是 广义 导数 应 用 发 展 中 具有 重要 意义 的 
HER, Heaviside 在 对 传输 线 理论 中 的 研究 中 引进 了 分 数 微分 的 概念 , 这 个 概念 由 Gemant 
所 推广 , 用 于 弹性 理论 的 研究 , 尽管 Heaviside 似乎 看 不 起 “ wet blankets of rigorists " ( JE 
谨 化 的 意义 )， 但 是 一 些 理论 家 意识 到 了 Heaviside 工作 的 优秀 方面 ， 并 且 试 图 在 为 人 所 
接受 的 数学 标准 方面 对 其 进行 改进 ， 1982 45, B.B.Mandelbrot 首次 指出 自然 界 和 许多 技 
术科 学 中 存在 大 量 分 数 维 的 事实 5-81, 并 在 整体 与 部 分 之 间 存 在 自 相似 现象 之 后 , 分 
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数 微 积分 成 为 了 研究 分 形 几 何 和 分 数 动力 学 的 有 力 工具 881-491, 并且 在 松弛 144-146), 
振荡 4 、 扩 散 和 输 运 理论 d9-5H 、 生 物 组 织 562-155| 、 高 分 子 材料 的 解 链 1521-96] | 
混沌 与 满 流 57158 、 随 机 游 走 99-60 、 统 计 与 随机 过 程 (2 -063 、 粘 弹性 力学 以 及 非 
牛顿 流体 力学 1971641165) 等 诸多 领域 得 以 应 用 ， 而 这 些 领域 的 应 用 反 过 来 又 促进 了 分 数 
微 积分 的 理论 研究 的 进一步 发 展 ， 成 为 了 当前 国际 上 的 一 个 热点 研究 课题. 

众所周知 , 微分 形式 在 物理 . 微分 几何 . 应 用 数学 领域 得 到 广泛 的 应 用 . 许多 偏 微分 方 
程 能 用 微分 形式 来 表达 . 近年 来 ， 通 过 将 其 建立 在 不 同 的 分 次 代数 的 基础 上 ， 在 文献 PA 
H, 外 微 积 分 又 被 拓展 另 一 种 企图 拓展 是 基于 非 结合 几何 而 建立 的 60, 最 近 ， Kathleen 
和 Mark 在 分 数 微 积 分 方面 作 了 大 量 出 色 的 工作 ，2001 年 ，Kathleen 和 Mark 给 出 了 分 数 
外 导数 的 定义 (67 ， 且 发 现 了 分 数 微分 形式 空间 为 有 限 维和 无 穷 维 向 量 空间 。 闭 形式 和 
恰当 形式 被 推广 到 新 的 分 数 形 式 空间 ， 并 且 在 特殊 情形 下 得 出 了 闭 的 可 积 条 件 . 此 外 还 
得 到 了 坐标 变换 规律 ， 陈 驶 教授 之 后 又 给 出 了 从 卡 氏 坐标 到 曲面 坐标 变换 2O. Kathleen 
和 Mark 又 对 分 数 微 积分 形式 作 了 进一步 的 研究 ， 特 别 地 考虑 了 分 数 微分 形式 空间 的 有 
限 维 子 空间 . 并 且 定 义 了 内 积 ， 偶 积 和 协 变 导 数 . 计算 了 积分 阶 形式 的 转换 规律 ， 并 用 矩 
阵 阶 分 数 微 积分 定义 了 抵 阵 阶 形式 从 而 产生 了 和 拢 阵 阶 形式 的 坐标 转换 规则 和 协 变 导 数 ， 
同时 还 指出 ， 除 了 矩阵 阶 是 非 对 角 情形 外 ， Poincaré 定理 对 所 有 阶 数 的 外 分 数 微 积分 都 
是 正确 的 . 


$1.3 ” 狐 立 子 理论 研究 概述 


孤立 子 (Soliton) 现象 最 初 是 由 英国 科学 家 Scott Russell 发 现 的 。 1844 年 ， Russell 
在 一 篇 题 为 《 论 波动 》 的 报告 中 记述 了 他 1834 年 观察 到 的 一 种 奇特 的 水 波 现象 . 当时 他 
正在 观察 由 两 号 马 拉 着 的 船 在 一 条 狭窄 的 运河 中 行驶 ， 船 突然 停止 了 前 进 ， 但 运河 中 被 
船 推 动 的 水 却 并 没有 停止 ， 而 以 测 涌 杰 腾 的 状态 聚集 在 船 头 ， 然 后 以 巨大 的 速度 滚滚 向 
前 ， 且 保持 着 巨大 的 轮廓 分 明 的 光 顺 孤立 的 妖 状 外 形 ， 显然 ， 它 不 改变 形状 与 速度 ， 沿 
运河 继续 前 进 。 他 骑 着 马 跟 踪 了 一 至 商 英 里 ， 在 运河 的 拐 咨 处 ， 这 种 孤立 行进 的 水 峰 才 
终于 消失 。 Russell 认识 到 这 种 水 波 现象 是 具有 关键 性 质 的 新 现象 、 新 事物 ， 随 后 进行 了 
更 加 细致 的 研究 ， 在 实验 室 作 了 很 多 实验 ， 用 多 种 方法 激发 ， 也 观察 到 了 同样 的 现象 . 
他 称 这 种 波 为 孤立 波 (Solitary wave) .但 限于 当时 的 数学 理论 各 科学 水 平 , 人 们 无 法 从 理 
论 上 给 予 这 种 现象 一 个 圆满 的 解释 ， 科 学 家 们 甚至 怀疑 孤立 波 现象 是 否 真 正 存在 ， 

在 随后 的 几 年 ， Airy Boussinesq 和 Rayleigh 等 人 相继 对 孤立 波 进行 了 研究 ， Airy 
得 出 结论 ， Russell 所 提 到 的 孤立 波 根 本 不 存在 ， Stokes 使 用 了 正确 的 方程 ， 却 得 到 了 错 
WER; Boussinesq 和 Rayleigh 分 别 从 数学 角度 证 明了 孤立 波 的 存在 性 。 Boussinesq 
为 近似 描述 孤立 波 , 提出 了 一 个 非 线 性 发 展 方程 ， 后 来 被 命名 为 Boussinesq FE. 但 是 ， 
Boussinesq 和 Rayleigh 的 工作 仍然 没有 使 那些 对 孤立 波 感 兴趣 的 科学 家 们 完全 信服 . 这 
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也 促使 荷兰 数学 家 Korteweg 和 他 的 博士 生 de Vries 对 水 波 现象 作 进一步 研究 . 

1895 Æ, Korteweg 和 de Vries 根据 流体 力学 研究 了 浅水 波 的 运动 ， 在 长 波 近似 和 小 
振幅 的 假定 下 ， 求 得 了 单 向 运动 的 浅水 波 运动 方程 ， 即 著名 的 Kav 方程 通过 一 定 的 数 
学 变换 ， KdV 方程 变 为 如 下 形式 


Uz + Guts + tsrs = 0 (1.3.1) 


其 中 % 为 波形 函数 .。 Korteweg 和 de Vries 从 上 式 求 出 了 与 Russell 所 发 现 的 孤立 波 现象 
一 致 的 结果 ， 即 具有 不 变形 状 的 脉冲 状 孤 立波 解 。 KdV 方程 的 解 ， 准 确 地 描述 了 浅水 波 
的 非 线性 特性 ， 行 波 速 度 依赖 于 其 本 身 的 振幅 ， 当 两 个 这 样 的 脉冲 波 沿 着 同一 方向 运动 
时 ， 波 峰 高 的 脉冲 波 的 行进 速度 快 ， 因 此 会 赴 上 前 面 波峰 低 的 波 而 发 生 碰 接 . 

然而 这 样 的 孤立 波 是 否 稳 定 ， 两 个 这 样 的 孤立 波 碰 擅 后 是 否 形变 ， 这 一 直 是 科学 家 
们 感 兴趣 而 又 无 法 证 实 的 问题 . 因此 在 没有 新 的 发 现 之 前 ，KdV 方程 以 及 孤立 波 仍 长 期 
处 于 被 埋没 之 中 . 

1965 年 美国 普林斯顿 ( Princeton ) 大 学 的 两 位 应 用 数学 家 M. D.Kruskal 和 N.Zabusky 
通过 数学 模拟 方法 深入 地 研究 了 等 离子 体 中 孤立 波 碰 擅 的 非 线 性 相互 作用 过 程 ， 他 们 总 
外 地 发 现 ， 两 个 孤立 波 在 磁 挤 后 各 自 的 波形 与 行进 速度 居然 都 能 保持 不 变 ， 仅 仅 是 相位 
发 生 了 玻 变 .这 一 性 质 使 人 们 联想 起 质点 粒子 和 波 粒 二 和 象 性 等 熟悉 的 现象 ,只 有 粒子 的 
碰撞 才 会 有 类 似 的 情形 出 现 ， 于 是 将 这 种 波 定名 为 孤立 子 Soliton )， 以 反映 其 粒子 属 
性 , 

“孤立 子 ” 没 有 明确 的 定义 ， 但 是 它 可 用 来 描述 一 个 非 线性 方程 或 非 线 性 体系 的 任 
意 解 ， 若 此 解 满足 ， 1. 可 表示 成 一 个 固定 形式 的 波 ; 2 是 局 部 的 、 误 变 的 或 在 无 穷 大 
时 变 为 常数 ; 3. 可 与 其 他 的 孤子 进行 强烈 的 相互 作用 ， 在 相互 作用 后 即使 于 加 原理 成 立 
其 形式 亦 不 会 改变 ， 

总 之 ， Kruskal 和 Zabusky 的 这 一 研究 工作 为 推动 孤立 子 理论 的 发 展 ， 树 立 了 一 个 
重要 的 里 程 碑 ， 此 后 ， 科 学 家 们 对 孤立 子 的 研究 兴趣 和 热情 便 一 发 难 收 ， 在 很 多 学 科 领 
域 都 发 现 了 孤立 子 运动 形态 ， 相 应 地 ， 在 数学 上 ， 发 现 了 一 大 批 共 有 孤立 子 解 的 非 线性 
发 展 方程 ， 而 且 已 逐渐 建立 起 较 系统 的 研究 孤立 子 的 数学 物理 方法 (69.707, 721 

目前 ， 较 为 完整 的 孤立 子 理 论 体系 正在 逐步 形成 ， 国 内 外 在 这 方面 出 版 了 很 多 专著 
(69, 71,72) 。 孤 立 子 理论 已 经 被 应 用 于 解决 等 离子 体 物理 、 凝 聚 态 物理 、 生 物 学 和 非 线 性 
光学 等 领域 中 某 些 难以 解决 的 问题 , 以 及 非 线性 作用 下 的 运动 规律 等 . 从 数学 方面 来 看 ， 
已 经 发 现 一 大 类 非 线性 发 展 方程 具有 孤 波 解 ， 求 解 方 法 也 出 现 了 许多 独特 的 分 支 ， 
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第 二 章 “分数 微 分 方程 以 及 分 数 微分 形式 
$21 “分数 微分 方程 及 其 应 用 


分 数 微分 方程 简介 
我 们 在 绪论 中 已 经 介绍 分 数 微 积分 是 指 任 意 阶 的 积分 和 微分 的 理论 ， 它 推广 并 统一 
了 整数 阶 微分 和 积分 的 概念 ， 让 我 们 来 考虑 n 重 积分 和 导数 的 序列 ; 


t t 
3 anf fnan, f fnan, 50, ŽO, fH. (2.1.1) 


任意 o 阶 的 导数 可 以 考虑 从 这 样 的 算 子 序列 引申 出 来 ， 并 由 Davis? 建议 记 作 ， 


D? f(8. (2.1.2) 


这 个 任意 阶 导数 也 简称 分 数 导 数 , 其 中 下 角 标 a 和 + 表示 分 数 微分 运算 运算 的 两 个 极限 ; 
也 被 Ross 称 作 分 数 微分 的 端点 69 . 同 祥 地 分 数 积分 表示 任意 阶 的 积分 ， 对 应 一 个 负数 
的 阶 数 ， 可 以 记 作 : 

aDË f(t), B > 0. (2.1.3) 


分 数 微分 方程 就 是 带 有 分 数 导数 的 方程 ， 分 数 积分 方程 就 是 带 有 分 数 积 分 的 方程 ， 
一 个 分 数 阶 的 系统 是 指 由 分 数 微分 或 分 数 积分 方程 描述 的 系统 【? 委 .在 分 数 微 积分 的 发 
展 史 中 , 关于 分 数 导数 和 积分 有 不 同 的 定义 , EBA: Griinwald-Letnikov 定义 ，Riemann- 
Liouville 定义 4 ， Caputo 意义 下 的 定义 [本 和 广义 函数 法 (76leon6) 。 限 于 篇 幅 不 将 一 
一 列举 . 为 了 介绍 分 数 数 积 分 和 微分 方程 ， 我 要 重点 介绍 Caputo 的 定义 . 虽然 Riemann- 
Liouville 定义 在 分 数 微 积 分 的 发 展 历史 中 起 到 了 相当 大 的 作用 ， 尤 其 是 在 纯 数 方面 ， 但 
是 它 逐 渐 不 能 满足 现代 的 科学 技术 要 求 ， 例 如 在 粘 弹 性 力学 和 固体 力学 中 ， 分 数 导数 是 
用 来 描述 材料 的 性 质 的， 而 且 建 立 在 流 变 学 的 基础 上 的 数学 模型 很 自然 地 引出 带 有 初 值 
条 件 的 分 数 微 分 方程 ， 这 些 应 用 的 问题 要 求 分 数 导数 的 定义 允许 使 用 带 有 实际 意义 的 初 
值 条 件 ， 例 如 未 知 函数 及 其 整数 阶 导数 的 初 值 条 件 ， 然而， Riemann-Liouville 导数 的 定 
义 下 , 未 知 函 数 分 数 导数 在 初 值 点 的 值 (或 极限 值 ) 是 常数 . 而 Caputo 导数 表示 的 分 数 微 
分 或 积分 方程 可 以 得 到 与 整数 阶 的 微分 和 积分 方程 相似 的 初 值 条 件 ， 也 就 是 说 未 知 函 数 
的 分 数 导 数 包含 未 知 冰 数 及 其 整数 阶 导数 (在 后 面 的 分 数 扩 散 - 波动 方程 中 可 以 很 好 的 
体现 ) 

Caputo 的 关于 分 数 导数 的 定义 [5 如 下 : 


g* f(t 1 t (n) 
Ee "rcg pott (n-1«Re(a)SmneN) — (214) 
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其 中 参数 a 是 求 导 的 阶 数 ， 可 以 是 实数 阶 其 至 是 复数 阶 ， a RBH f 的 初 值 {在 本 文中 
将 只 考虑 a 为 正 实数 阶 的 情形 )， 在 Caputo 导数 的 定义 下 的 我 们 可 以 得 到 ; 
TE 
ate 
这 一 点 不 同 于 Remann-Liouville 的 定义 ， 因 为 在 R-L 定义 下 的 分 数 导 数 对 常数 求 导 往往 
不 等 于 0 . 


= 0, (C is a constant) (2.1.5) 


B "T 
a = 2 Mosen (2.1.6) 
öt tf-0, (852-1) 


FELD 
类 似 于 整数 阶 导数 的 性 质 , 当 f(r) 在 [ot] 连续 , H (7) 在 [atl A n+1 阶 连续 导数 时 ， 
Caputo 分 数 导 数 同样 满足 线性 和 Leibniz 法 则 ; 


ag AQ) Holt) = Ages (t) + pao (11) 
(k) ge-k 
Z enro) = a Je. (2.18) 
其 中 RH. 
Caputo 分 数 导数 的 Laplace 变换 和 Fourier 变换 为 : 
(Ds) gut - Yr eto, (219) 
k=0 
和 " 
Fl aah (05) = (Ci) (0 (0i). (2.1.10) 


求解 一 般 的 分 数 微分 方程 的 一 种 最 有 效 的 方法 是 Laplace 变换 法 ， 它 是 基于 Laplace 
变换 技巧 和 双 参 数 的 Mittag-Leffler 函数 (参见 下 一 节 ) 的 Laplace 变换 公式 的 一 种 方法 . 

Mittag-Leffler 型 函数 

我 们 知道 指数 函数 ez 在 整数 阶 的 微分 和 积分 方程 中 起 到 了 十 分 重要 的 作用 . 同样 ， 
它 的 参数 推广 ， Mittag-Leffler 型 函数 在 解 分 数 微分 和 积分 方程 中 也 起 到 了 重要 的 作用 . 
M-L 型 岗 数 很 自然 地 起 源 于 积分 方程 的 解 ， 尤 其 是 广泛 地 应 用 在 动力 方程 ， 随 即 游 走 ， 
levy 飞行 , 和 所 谓 的 超 扩 散 运 输 方程 的 分 数 推广 方程 中 . 近来 ，M-L 型 函数 的 重要 性 逐渐 
被 认识 , 而 且 被 更 多 地 关注 . 不 仅 是 从 解析 的 观点 ， 而 且 从 数值 分 析 的 观点 ， 还 包括 用 于 
描述 分 数控 制 系统 ， 分 数 粘 弹性 模型 等 等 很 多 方面 .它们 的 定义 以 及 性 质 在 最 近 的 关于 
分 数 微 积 分 , 分 数 扩散 方程 和 积分 方程 , 力学 方面 的 书籍 和 综述 中 可 以 找到 88. 74 7) | 

Mittag-Leffler 是 一 个 整 函数 ， 在 T 中 定义 如 下 : 


E,(z) = ea >0 (2.1.11) 
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E e =F rim 的 单 参数 推广 ,同样 地 ， MI 再 函数 的 双 参 数 的 定义 T y 


Bglz) = YE {a > 0,8» 0) (2.1.12) 
自然 地 我 们 得 到 了 如 下 结果 : 
Bays d ye oe 2.1.13 
u= d rj AET . — (2143) 
Ei2(z) = » 2 3 a = PUT (2.1.14) 
£.T(*2) Skt z 

BaG)-Ymies-) p m (2.1.15) 

PR +3) (k +2)! z2 o ' a 
f T P-E plate)dt =. (el< 1) (2.1.16) 

0 l-z 


令 8=1 我 们 可 以 得 到 函数 的 单 参数 形式 : 
Eji(z)- V = = EB,(z). 2.1.17 
1 By ak 7 (z) (2.1.17) 
为 了 求解 有 理 阶 微分 方程 ， 文 [6 提出 了 如 下 函数 


e, — t" » T mns = P E v41 (at). (2.1.18) 

Mittag-Leffler E83 £4) 25064 75 $2 PHA 

我 们 知道 扩散 方程 ， 波 动 方程 是 数学 物理 方程 的 基本 方程 ， 它 在 多 门 学 科 中 都 有 重 
要 的 应 用 . 随 着 应 用 的 需要 ， 在 物理 ， 金 融 ， 光 学 以 及 生物 学 等 诸多 方面 出 现 了 用 带 有 
任意 阶 时 间 和 空间 导数 的 分 数 微分 方程 描述 的 模型 . 其 中 最 重要 的 基于 Lévy 随机 游 走 和 
分 数 布朗 运动 的 分 数 扩散 方程 ， 虽然 对 于 一 系列 的 分 数 扩散 - 波动 方程 (FDWE) 已 经 有 
了 十 分 合理 的 物理 解释 ， 但 是 它们 的 初 值 和 边界 向 题 的 解 并 没有 给 出 .为 了 解决 这 个 问 
题 - 找到 所 谓 的 基本 解 ， 我 们 在 高 维 空间 中 定义 了 推广 的 ML RHA CB 
多 的 重要 的 积分 变换 、 同 时 得 到 了 一 些 有 用 的 结果 . 

这 种 新 的 M-L SURGE MOF: 


e(t. y; o, 8,7) = 197 Eag- D|yl't*), (2.1.19) 


其 中 + 代表 时 间 变 量 ， y 代表 空间 变量 ， 记 作 : y= (yy yn) € R BY 可 以 为 任 
意 实 数 阶 且 D 是 有 一 定 物理 解释 的 物理 常数 . 
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我 们 发 现 这 种 形式 的 函数 已 经 频繁 地 应 用 在 求解 动力 系统 ， 分 数控 制 ， 和 分 数 扩散 - 
波动 方程 的 Laplace 方法 中 . 然而 由 于 没有 现成 的 公式 表 和 手册 可 以 查 到 该 函数 的 Laplace 
变换 ， Fourier 变换 以 及 相关 变换 ， 而 且 计算 非常 复杂 ， 我 们 给 出 了 它 的 一 些 重要 的 积分 
变换 ， 

函数 c(t, yio, 6,7) 关于 时 间 变量 + 的 Laplace 变换 为 

setti Bip) = f P-E gl- Dirty 
[9 DIMT y 
a purs EE (ak 4- 1) 


Den f ent yak At 


-D 
Fa (ok Bp 65 
p8 


= p+ Djy 


e(t, yia B, y) 函数 的 关于 时 间 变 量 t 的 Mellin 变换 可 以 通过 公式 


A (f)s) = rü- mig pU) jphl- s) 


得 到 
A (e(t yia. B, y) s) = M CE (e(t uio 8 Y)ipy;l- s) 


pi. dos 
NI DBP 


(2.1.21) 
ra p=) / y : spp 


_ (Db) B+s-1, , 1-8-5 
—. al(1- s) i qo mg ) 


其 中 BAS 1+ Es) 是 Beta BR. 
e(t, y; o, B, y) 函数 的 关于 高 维 空间 变量 的 道 Fourier 变换 十 分 复杂 ， 


E(t, sa, 8,7) = F Helt, ya, p, y) z} = dus 人 . € *Ve(L y; a, B, y)d" y, (2.1.22) 


其 中 是 Ty 高 维 空间 变量 ， infe: t= (21,22, enn) V md (91, Y2, Yn) cR" 
上 式 很 难 用 通常 的 积分 变换 直接 得 到 结果 .所 以 我 们 利用 Melin 变换 和 特殊 函数 的 技巧 
对 它 进行 约 化 
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关于 时 间 上 ， 我 们 先 对 上 式 的 两 边 进行 Mellin 变换 


Miel 2:0, 0, y) s) = FTHM lelt, zia 8,7); 8}; 2} 


-1 (Diy M -1 1-8- 
i A s) BÊ 7 Me La) (2.1.23) 


BUR psc). Ces 


-aaa ht 9 P f) 


其 中 * 是 一 复 变量 , 满足 51 < Re(s) < s2 (51,52 是 实数 ) .对 空间 变量 y = (y1, Y2 Yn) € R° 
作 正 交 变 换 和 球 坐 标 变换 ， 我 们 可 以 得 到 : 


Qn) 2 ["? a, Ip-s 
z} = es | peo Ja_1(p|x1) 40s (2.1.24) 


aü-£-2 


$^ 


其 中 p = [yl Jg- 是 第 一 类 (3 — 1) Bt Bessel 函数 .利用 第 一 类 Bessel MAHA 


1 ONE) 


质 J v1, bdt = s EE 我 们 可 以 得 到 ， 
Jaspi 225 a0-6-5) 
$7 (w x ua C MEME De fer ra ER (2.1.25) 


将 (2.1.22) 式 代入 (2.1.20) 式 我 们 得 到 : 
(PD) p(n + erg + 1872-2) 


ME (t2; a, B, y); 8} = ————ÀHu 2.1.26 
{el Tj B y) s) an? sp T s) rü — s)r( nere- 22) ( ) 


FEM ERGIT Mellin 变换 产生 ， 


e(t,2;0, 8,7) = £^ (e(t. y; a, B, y); 2} 
ol De pete fa (At hr + i-r) (n + 2078-2) 
"VD" E Bri iw 2D rü- gy EE) 
其 中 c < Re(c) < ole ceo 是 实数 ); 上 式 中 的 积分 路 线 是 Bromwich R44. WH, Fox's H 
A ( 详 见 781 的 附录 ) 的 定义 ， 可 以 知道 : 


一 — $71 (DU eol zy 0-822203) 
£(5,210,8, y) = at yma) BP fair) ^ Psp lo- AE UD) 0-H, p 
(2.1.28) 


其 中 H (ale Ic ets hg gen xoc az, g,) 是 Fors 本数 ,根据 文人 附录 中 的 
公式 (A.7)(A.8) 可 以 写作 : 


(2.1.27) 


x i971 zl" 
£(t,z;a, B, y) = S (e(t, y; a, B, vy); z)p- = NS uc. i ETENI (2.1.28) 


下 面 看 看 新 构造 的 这 一 类 函数 的 应 用 
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通常 的 分 数 扩散 - 波动 方程 可 以 表示 为 : 


Ou 
ET 


其 中 A3 ARIK Laplace 算 子 ， 在 (80,20) 中 第 一 次 利用 Fourier 变换 给 出 定义 ， 


--D(- A)?u (0<a,7 x 2), (2.1.30) 


F{(-A) ult, 2; o, 7);w} = wl (u(t 2; a, y);w} 


其 中 D 代表 一 个 物理 常数 ， ,8 可 以 是 任意 实数 ， 相 关 的 工作 可 以 参见 (81, 82 83,84, 

现在 我 们 在 n 维 空间 (n > 1) 内 考虑 初 值 - 边界 问题 ， 

sult tia, y) = ~D( —A)?u(t, z; a. y), (£202z€cR'"m-i«c«aetm) (2.1.31) 
T GUN u(t, x; o, y) = 0; o; Fue zi2T) = felz). (k 20,1,...,m - 1) (2.1.32) 
其 中 的 分 数 时 间 和 空间 导数 都 是 Caputo 定义 下 的 导数 . 

考虑 边界 条 件 ， 方 程 两 边 同 时 对 时 间 变 量 t 进行 Laplace 变换 ， 对 空间 变量 z 进行 
Fourier 变换 可 以 得 到 ， 

b» Fy(u)po-1-* 


U(p,u; a, Y) = 


其 中 U (p, wa, y) 是 ult, sia, y) 的 Fourier-Laplace 变换 。 Fy (w) 是 fic); p 的 Fourier 变 
i&, pÆ Laplace 变换 参数 ， w 是 Fourier 变换 参数 .利用 (2.1.20) A, AJAH: 


m-i 
u(t, wa, y) = 35 Flw)elt,w; a, k 4-1, 9). (2.1.34) 
k=0 
其 中 u(tu;a, y) 是 u(fz; o. 7) 关于 参数 w 的 Fourier 变换 。 再 对 上 式 进行 道 Fourier BH 
可 以 得 到 方程 (2.1.31)-(2.1.32) 的 解 : 


m=i 
uena = $5 | AOGE- (2.1.35) 
k=0 R” 


其 中 
Gí(t,z) = ^7! (e(t, a, k + 1,y);x). (2.1.36) 
正如 整数 阶 的 情形 ，Gi(t,z) 也 可 以 叫做 分 数 扩散 - 波动 方程 的 基本 解 , 利用 (2.1.29) 
式 我 们 可 以 得 到 ， 


ean (he lm i 


Gi (t, 7; o, y) = = Us Hs (aye (1,1),(§ 53). F)”? 
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进一步 地 ， 利 用 特殊 函数 的 人 性质 我 们 可 以 得 到 数学 物理 方程 的 很 多 经 典 结 果 . 
1422: (2.131)-(2.1.32) 可 以 约 化 为 : 


Su za, y) = DAu(t, z; o. 7), (t»0zeR'm-l«aoetzm) (2.1.38) 
k 

lim u(t,z;o,y) = 0; Zua z;a,'y) = fu(z). (k 2 0,1,.., m — 1) (2.1.39) 
jz|-^--oo ot 


这 就 是 带 有 分 数 时 间 导 数 的 分 数 扩散 - 波动 方程 55] ， 直 接 从 式 (2.1.37) ， 我 们 可 以 可 到 


ew tala EP eie) (EO Lum 1) (2.140) 


Gk (t, 250,2) = UR ADE 04.1) 


这 个 结果 与 文 991 中 的 结果 一 致 . 
特别 地 ， 当 n = 1 时 ， 我 们 得 到 ; 


EFE W(-- 坟 人 j 1+k- 2). (k 20,1,.., m — 1) (2.1.41) 


它 与 文 189147] 研究 得 到 的 结果 一 致 
k=0 , RATA: 


Gi(t, 2; a,2) = 


EL. mu (2.1.42) 


1 
Gi 1,2;0,2 三 a en] , 
ol = BPS r2 


它 可 以 被 称 作 分 数 Green 函数 19 . 
当 a= 1 时， 上 式 就 可 以 约 化 成 经 典 的 表达 式 : 


1 dz 


E je]? 
Go(t,2;1,2) = 2J/Di -Fi -55 ts Ue zp(— ib) (2.1.43) 
它 就 是 我 熟知 的 扩散 方程 : 
ERI = DAu(t,2;0,7). (2.1.44) 


的 Gaussian 分 布 解 ， 参 见 (891. 
当 a=28= 1 时 ， 丛 (2.137) 我 们 又 可 以 得 到 经 典 的 波动 方程 的 基本 解 


e 0, if n is an odd greater than one, 
t,2;2,2) = i 2.1.45 
zc ame n... if n = 2l is an even where l € N ) 
vD(21-1)21) (Dt? -|z) 2 
21. 方程 (2.1.31)-(2.1.32) 又 可 以 化 简 为 : 
8 
apum) -D(-AfMwuwtz;y,  (t»0zeR) (2.1.46) 
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k 
lim u(t,z;1,7) =0; 2 ulg z;1ly)- felz). (k=0) AF, (2.1.47) 
|z| 一 十 ce at 


该 方程 是 通过 连续 随即 游 走 模型 构造 的 ， 用 来 描述 混沌 和 汪 流 中 的 超 扩散 现象 的 方程 
进一步 考虑 简单 的 情形 7 = 1 = 0 ， 利 用 Gamma RRA FoxH 冰 数 的 性 质 我 们 可 以 得 
到 波动 方程 的 Gauchy 分 布 解 : 

Golt, x;1,1) = : IT( 吐 )vZ 


RM ee 2.148 
nF (|x? + D2:22) F n. 


它 与 文 (86,871 中 的 结果 相符 . 
值得 指出 的 是 ， 当 0 < a < 1 是 方程 (2.1.38) 表示 所 谓 的 慢 扩 散 ， 当 1 < a < 2 时 被 
称 作 中 间 过 程 ， 这 些 方程 都 有 十 分 明确 的 重要 的 物理 意义 . 


82.0. 分数 微分 形式 及 其 应 用 


微分 形式 

微分 形式 的 微 积分 是 纯 数 学 的 一 个 很 讲究 的 分 支 , 是 应 用 数学 的 一 个 强 有 力 的 工具 . 
关于 强调 应 用 方面 的 领域 ，Flandersl83i 给 出 了 简明 的 介绍 . 下 面 给 出 微分 形式 的 一 个 简 
短 的 回顾 .向 量 空间 在 一 点 pe E" (on 维 欧 几 里 德 空间 ) 能 被 构造 出 下 面 的 表达 类 型 : 

一 形式 ， | 


a= Y adz; (2.2.49) 
i=l 
二 形式 ， . 
B= 2 bi;dzi ^ dz; (2.2.50) 
$jzi 
“n ”形式 ， 
w = wdr1 ^ dr2... A dEn (2.2.51) 
其 中 {zi} 是 B? 的 卡 氏 坐标 ， 上 面 的 和 式 是 对 满足 条 件 
dzi A dzj 一 —dz; A dzj (2.2.52) 


的 所 有 指标 进行 的 函数 a 和 bs 等 仅 依赖 于 p ， 并 且 依据 不 同 的 应 用 可 以 取 实 数值 和 
复数 值 . 如果 一 个 大 形式 7 乘 一 个 mm 形式 x ， 则 下 面 是 正确 的 : 


y^g — (C1) uay (2.2.53) 


AUR k+m>n, 结果 是 0 。 外 积 是 分 配 、 结 合 和 反对 称 的 .容易 看 出 ， 有 形式 的 向 量 空 
闻 的 维 数 是 ; 
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n n! 
(i) m 4 k >n BH 0. 


通常 的 外 导数 的 定义 为 ] 
d- Y) duc. (2.2.54) 
i=l li 
外 导数 将 形式 映 入 上 十 1 形式 ， 且 有 下 面 的 代数 性 质 . 让 7 和 入 是 上 形式 ，/ 是 一 和 nm 
形式 ， 那 么 


d(y + A) = d(y) + d(4) (2.2.55) 
d(y Ap) = (dy) A p+ (-1) y ^ du (2.2.56) 
d(dy) = 0 (2.2.57) 


最 后 一 个 恒等式 被 称 为 Poincaré 引 理 . 形式 y 称 为 闭 的 , WR dy = 0 . 形式 y RAK 
当 的 ， 如果 存在 形式 u 满足 du =q, AP u RE y 的 阶 数 小 1, 恰当 形式 一 定 是 闭 的 ， 
但 闲 形 式 不 一 定 是 恰当 的 ， 有 兴趣 的 读者 可 以 参阅 Flanders? 或 者 Lovelock 和 Rundl89 


得 到 更 详细 的 叙述 . 
分 数 微分 形式 空间 
如 果 在 外 导数 的 定义 中 导数 允许 取 分 数 阶 ， 那 么 就 能 定义 一 个 分 数 外 导数 107 
n A ov 
d'z Yo az} (8(z; - a) (2.2.58) 


i=] 


注意 下 标 i 表示 坐标 数 ， 上 标 v 表示 分 数 坐 标 微分 的 阶 数 ， ai 是 导数 的 初始 点 (在 
本 文中 考虑 a; = 0 的 情形 ) . 


有 时 候 用 标记 OP 表示 
(O(a; ~ a;))” 
在 二 维 倩 形 (zx,y) ， 取 初始 点 为 原点 ， 则 ze 的 v 阶 分 数 外 导数 为 
v T(pc 1) p-u v a? 
对 导数 参数 的 特殊 值 可 以 得 到 如 下 的 结果 ， 
v=0, dx? = 22? (2.2.60) 


v=], dlz?-dzlpz?! (2.2.61) 


15 


RE: 分数 微分 方程 和 非 线性 发 展 方程 中 问题 的 研究 


v—2, dg? = dz?p(p - la"? (2.2.62) 


和 标准 外 微 积分 类 似 ， 用 dry 可 以 构造 出 向 量 空间 . $ F(v,m,n) 是 一 个 在 pe E^ 
的 向 量 空间 ，，* 表示 基 元 素 的 分 数 微 分 阶 数 的 和 数 ， m 表示 出 现在 基 元 素 里 的 坐标 微分 
的 个 数 ，n 表示 坐标 个 数 ， 并 且 {zi} 是 pe E^ 的 卡 氏 坐标 , 例如， F(v1n) 的 一 个 基 
序列 是 {dx}, dzj,...,dzz) 并 且 Fo, Ln) 的 任意 元 素 可 表示 为 


a= Y aidz? (2.2.63) 
i=l 
NAEK v, 这 是 一 个 n 维 向 量 空间 .顺便 指出 ， 对 每 一 个 v 的 值 存在 一 个 不 同 的 
向 量 空间 . 当 v = 1 时 ， 重 新 得 到 来 自 外 微分 的 一 形式 . 现在 假设 基 元 泰 由 两 个 坐标 微分 
构成 ， 即 ftv,2,n) ， 在 这 种 情况 ， 基 序列 是 很 复杂 的 ， 


{act Adak? [ij € (1,2,..., n], Vus pa 2 0, p1 + ye = v} (2.2.64) 


注意 到 dz Ade}? 是 0 SHARA m =m, Hoi j 等 等 , 故 ，F(v,2,n) 的 任意 一 个 元 素 
可 被 表示 成 以 下 形式 的 和 : 


B= xx (Bij (vis 了 一 vi)dz; dr; ")dwi (2.2.65) 
i=1 j2179 ; 

其 中 fun) =0。 和 前 面向 量 空间 P) 不 同 ， Ptw 2.) 对 任意 的 v 的 人 都 是 无 穷 
维 的 ， 而 且 是 不 可 数 无 穷 维 的 . 进一步 ， 类 似 于 标准 的 外 微分 的 情形 ， 我 们 可 以 构造 空 
[B]: 

F(v,m,n) = {de pda... A dat [pg m v) (2.2.66) 
Rt 表示 基 元 素 的 微分 阶 数 且 可 以 为 任意 阶 ，“m;' 表示 坐标 数 ， 空 间 中 的 元 素 就 叫 
fit v- BA. 容易 看 出 Flv,m,n) 是 不 可 数 无 穷 维 空间 当 m > 2 时 , $ A € Fm, n), 
BeF(ukn), WAR 4 和 B 具有 反对 称 性 质 ， 即 : 


AAB= (-A"BAA€ F(utu,k+m,n) {2.2.67) 
不 同 于 标准 外 微分 形式 的 外 代数 ， A 人 8B AERP, Á ktm>nij, d $Æ F(u k,n) 


PERRA FA 十 vk 十 1,n) m. 
为 了 研究 F(vm,n) 的 代数 性 质 ， 先 来 考虑 它 的 有 限 维 子 空间 ， 


m 
G(m zz2 Vm, n) = [[em.n) (2.2.68) 


izl 
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其 中 G(y,n) = Fv,1,n) H dimG(vi,vo,--- um n)=n™., WY 


F(v,m,n) = U GÒ, ¥a,+++ Vn) 
T nsv 
i=l 


接 下 来 , 在 我 们 比较 熟悉 的 有 限 维 空间 Gi, ww,… ,vm,n) 中 , 我 们 可 以 构造 内 积 ，Hodge 
积 而 且 所 谓 的 分 数 Poincaré 引 理 也 成 立 : 


dd (a) —0,Va € F(v,m,n), WER. (2.2.69) 


更 详细 的 内 容 可 以 参见 68, 
我 们 知道 梯度 ， 旋 度 和 散 度 是 向 量 分 析 中 的 基本 求 导 算 子 ， 当 引入 了 外 微分 后 它们 
的 形式 可 以 得 到 统一 。 近来， 它们 也 被 推广 到 了 分 数 导数 的 形式 80, 91] ， 而 且 它 们 被 应 
用 在 解决 物理 问题 当中 ， 有 一 定 的 物理 解释 .基于 以 上 的 事实 我 们 考虑 用 分 数 外 导数 来 
引入 分 数 梯度 ， 旋 度 和 散 度 . 
WH RS : (x1, 22,23) EK KAA. 
1. RR f Æ R 中 的 光滑 函数 ， 按 照 定义 (2.2.58) 我 们 可 以 得 到 ， 
„e Of a of 
f= Gaye t BY (za) 
则 我 们 称 Gro ade de) 是 了 的 分 数 梯度 ， 记 作 gradj, 
2. Hw” = asda 十 a2dzy + agdz¥ € G(v,3), oi(i = 1,2,3) Æ Ra 上 得 光滑 函数 ， 则 利用 
(67) 中 的 乘积 法 则 (18) ， 我 们 可 以 得 到 ， 


dal, (2.2.70) 


dz 十 


3 
dw" = X "d" (adz) = Yu ^ AM Jt r Ba, =k mae aem. (2.2.71) 
i=l ij=1 
注意 到 
E 
ar (d) =0,Vk>1 
(2.2.71) 可 以 约 化 为 : 


e" 
Pw” = Sar ^ daz ie Gaye 2 ^ dz; Gey y ~ Ome! 


17 一 工 

vag — Om 03 iv jw (2.2.72) 
= loco is izle ^ del + y 7 (5s yl s Adz} 

à" ag ov a) 


aye (as (8s; 771 ^ das. 


我 们 称 (合生 — Goi. ed - gu Bor duke) 是 向 量 场 X (a1, 02,03) BNE 


度 ， 记 作 curl’X , 
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3. 设 必 =aliadraAdadrs+oazdzsAdzl +asdzl A dra € G(v,v,3) ， 因 为 


3 
d" (a1dz ^ dz) = > dr% ^ — (adu; ^ day) 


E 


-Me M eee (aya a e Ade 


c y 


而 且 注 意 到 : m 
gi nt ^Adz)-0, Vk>1, 


我 们 有 ， 


3 
: e Br 
j=l 


因此 有 


d -一 一 03)dz1 A dz ^ dz 


" » 
E a a 
HRR, RI Gap + Gea + hune 是 X = (a1,02,03) 向 量 场 的 散 度 ， 记 作 : div X 

在 情形 1 2, 3 中 , v 可 以 使 任意 实数 甚至 是 分 数 阶 的 , 当 v = 1 时 ， 所 谓 的 分 数 
梯度 ， 旋 度 和 散 度 就 可 以 约 化 成 经 典 的 向 量 场 中 标准 的 梯度 ， 族 度 和 散 度 . 此外， 类似 
于 Poincaré 引 理 ， 由 分 数 Poincaré 引 理 也 可 以 得 到 : 


curl"grad"f =0,  div"curl" X = 0, (2.2.73) 
当 v=1 了 时， 又 可 以 得 到 经 典 的 结果 : 
curl(gradf) =0, div(curlX) = 


值得 提出 的 是 , 我 们 提出 的 grad” f, curi" X, 和 div’ X. 与 文 P0. 9 中 用 其 他 试验 和 建 模 方 
法 得 到 的 结果 一 致 。 而 且 这 些 工具 将 进一步 推动 分 数 微 积分 的 发 展 ， 而 且 可 以 更 好 的 解 
释 许 多 物理 现象 ， 例 如 分 数 布朗 运动 和 随即 游 走 ， 


大 连理 工大 学 硕士 学 位 论文 
第 三 章 ”“ AC 一 BD ”理论 与 C-D 对 的 构造 方法 


自从 张 鸿 庆 教授 821 于 二 十 世纪 六 十 年 代 担 出 了 “ AC=BD ”思想 ， 并 于 1978 年 发 
表 以 来 ， 他 和 他 的 学 生 们 在 这 方面 做 了 大 量 的 工作 ， 使 得 这 一 思想 在 电动 力学 、 弹 性 力 
学 流体 力学 、 量 子 力学 、 孤 立 子 理论 、 物 理学 等 方面 得 到 了 广泛 的 应 用 . 这 一 思想 是 一 
TARA, BG NMA. RH. BA” 的 原则 .近年 来 该 思想 推广 到 解决 非 线 性 问题 
中 ， 张 鸿 庆 教 授 又 提出 了 OD 可 积 系 统 与 C-D 对 的 概念 ， 形 成 了 在 微分 方程 (组 ) 求解 
中 的 C-D 可 积 理 论 ， 在 孤立 子 理论 及 其 应 用 方面 有 了 很 好 的 成 绩 ， 本 章 简要 介绍 张 鸿 庆 
教授 提出 的 关于 微分 方程 (组 ) 求解 的 “AB=CD ”理论 及 应 用 ， C-D 对 的 构造 方法 


$3.1 “AC=BD ”理论 及 其 应 用 


“AC=BD "理论 的 基本 思想 就 是 将 复杂 不 易 求解 的 方程 ( 原 方程 ) 通过 适当 的 变换 

转换 为 简单 易于 求解 的 方程 ( 目标 方程 ). 不 失 一 般 性 ， 可 形式 地 将 原 方 程 和 目标 方程 分 
别 表 示 为 hv = 0 和 Dv =0， 则 原 方 程 的 求解 就 变 为 寻找 适当 的 变换 uw = Cv ， 将 原 方程 
化 为 易于 求解 的 目标 方程 Dv = 0 . 但 是 , 在 实践 中 往往 需要 求 得 算 子 B. ( 辅助 算 子 ) 使 
其 满足 “ AC=BD "， 有 时 还 需要 求 得 算 子 R ( 余 算 子 )， 使 得 AC-BDER. 其 具体 格 
AE: BR 4u = 0 为 待 求解 的 原 方程 ， Dv = 0 是 易 求解 的 旭 标 方程 ， 寻找 变换 = Cv 使 
f$ Au 20— Dv 20, H CkerD = ker4 .对 一 般 微分 方程 的 求解 ， 就 转化 为 以 下 问题 的 
解决 : 给 定 算 子 A ,构造 算 子 C 和 D ， 使 得 CkerD = ker4 ， 及 如 何 构造 变换 = Cv, 
将 待 求解 的 方程 Au = 0 约 化 为 目标 方程 Dv=0. 

定义 3.1.1: RX BREESE, A,B,C,D ÆA X EX 的 算 子 ， 对 任意 ve X， 

AC (v) = A (Cv), BD (v) = B (Dv) 

如 果 对 veX ，4cv = BDv， 则 称 4C = BD. 

定义 3.1.2: 如 果 对 于 算 子 4 , 存在 算 子 B,C, D , 使 得 4C = BD, CKerD = KerA, 
其 中 KerA = {u | Au = 0), KerD = {v | Dv = 0}, WR Au = 0 是 可 积 系统 . 
d CKerD# KerA, (B CKerD C KerA WK Au = 0 为 部 分 可 积 系 统 ， 


定义 3.13: 算 子 C 和 D 称 为 算 子 4 的 C-D 对 ， 如 果 系 统 : 


| Civ (3.1.1) 


D(v,u) = 0. 


其 相 容 条 件 恰 为 Au=0, EF u AAR, “ 恰 " 的 意义 为 如 果 系 统 ( 2.1 ) 的 另 一 个 兼容 
UE A*u m0, ABA KerA*c KerA, di C(v,u) =0 可 写 为 w=Cv，D({v,w) =0 可 写 


18 


RR: 分 数 微分 方程 和 非 线 性 发 展 方程 中 问题 的 研究 


为 Dv =0,， 那么 4 有 显 式 的 C-D 对 . 若 CKerDC KerA, JW Dv =0 HER v, 
Hu=Cv, Wl] Au- 0, 若 CKerD > KerA, 则 对 hwu=0 的 任意 解 u, 必 有 ve KerD, 
使 得 w= Cv. 如 果 CKerD D KerA fü CKerD C Kera 同时 成 立 ， 即 CKerD = KerA, 
这 时 方程 hu = 0 的 一 般 解 为 = Cv ， 其 中 v 满足 Dv = 0 ， 也 称 在 变换 = Cv 下, X 
程 Au =0 5 Dv =0 等 价 . 


定义 3.1.4: ABA foo : = h 的 相 容 性 条 件 为 Au 20, 则 称 hu =0 是 C-D 
U,uj = U. 
RH, FB one de. J3 Au =0 BY C - D 对 . 
D(v, u) = 0. 


定理 3.1.1: RX 是 线性 空间 ，A, B,C,D EX BX 的 线性 算 子 . 如 果 4C = BD， 
B(0 20, B CkerD > ker4 ， 则 方程 (组) Au = 0 的 一 般 解 为 4 = Cv, KH v W 
足 Dv=0. 


定理 3.1.2: 设 
01] Q12 ... Qin 
0231 4022 ... Qn 
A= 
Qn) ûn? ... Onn 


其 中 o 是 线性 偏 微分 算 子 ， B,C,D 是 偏 微分 算 子 矩 阵 ， 且 满足 AC=BD, CKerD= 
Ker4 ， 则 非 齐 次 方程 Au = f 的 一 般 解 可 表示 为 4 = Coe, Dv-g, Kt eg 是 方 
程 Ae + Bo = 的 一 组 解 . 

推论 3.1.1: 车 和 是 线性 空间 且 CKerD 2 KerA ， 则 Au = 0 的 解 可 以 用 us = Cv, 
ok, KH on WE De, =0. 

推论 3.12: RX 是 线性 空间 ， ABCDEX AX HARERT, fex, 
HAC=BD, CkerD 5 kerA ， 则 方程 (组 ) Au = f 的 一 般 解 为 ， w= Co +e， 
Du=9., 其 中 e.g 满足 方程 (组 ) Ae-Bg-f. 

证 明 : WREE XA X HAT M,N H EE AM -BN = E, lle = Mo I g = Nó 
满足 方程 he + Bg = f ， 其 中 满足 方程 B8 = 了. 

张 鸿 庆 教 授 及 其 学 生 应 用 “ AC=BD ”理论 在 微分 方程 求解 方面 作 了 大 量 的 工作 ， 
这 一 思想 是 开放 性 的 ， 所 以 在 很 多 方面 得 到 了 广泛 的 应 用 ,如 弹性 力学 、 电 动力 学 、 流体 
力学 、 量 子 力 学 、 孤 立 子 理论 、 理 论 物 理 等 . 目前 ， 算 子 A 和 DD 可 以 有 如 下 不 同 的 表达 
形式 : 


20 


大 过 理工 大 学 硕士 学 位 论文 


A — D 
任意 微分 方程 组 具有 对 角形 式 的 微分 方程 组 
非 线 性 微分 方程 线性 微分 方程 
变 系 数 微分 方程 常 系数 微分 方程 
高 阶 微分 方程 低 阶 微分 方程 
高 维 方 程 低 维 方程 
微分 方程 代数 方程 
任意 的 方程 具有 特定 性 状 的 方程 
不 可 分 离 变量 微分 方程 可 分 离 变 量 微 分 方程 
不 会 求解 方程 会 求解 方程 或 具有 重要 性 的 方程 


83.2. C-D 可 积 系统 及 其 构造 方法 


定义 3.2.1 R Au = 0 为 一 个 具有 解析 系数 的 微分 方程 系 ， 如 果 它 有 形式 寡 级 数 解 ， 
则 称 它 是 形式 可 积 的 ， 如 果 这 个 级 数 收 剑 ， 则 该 解 是 古典 解 . 

定义 3.2.2: BA RATER BAER EAR EAR BOR A 
况 下 可 积 条 件 的 总 称 . 


定义 32.3 车 方程 组 n Vey 


D(u,v) = 0. 


的 相 容 性 条 件 为 Au = 0, 则 称 Au = 0 是 CD 可 


s. XH s = 0, 
D(u, v) = 0. 

最 近 Reid 提出 了 将 一 类 微分 方程 约 化 为 标准 型 ， 进 而 确定 解 空间 的 维 数 ， 最 后 给 出 
T Taylor 形式 震级 数 解 的 思想 和 算法 (详细 内 容 请 参看 文 99), 

对 于 一 个 含有 m 个 自 变量 x = (zzm) 及 n THEA VH), V^(x), V? 的 导 
数 表示 为 DoV? = 3%t*…tomV?/8z9.…68zom, FH a = (az, yam) € NmN = {012 
ord(a) = a 十 oa 十 … 十 am 之 0 叫做 导数 的 阶 . | 

首先 定义 一 个 序 ， 我 们 说 Dav? > DV, 如 果 满 足下 面 的 条 件 ， 

(1) ord(a) > ord(b); 

(ii) ord(a) = ord(b) H p >q; 

(ii) ord(a) = ord(t),p = q, E a 8 的 第 一 个 非 零 分 量 为 下 

在 上 述 序 下 ，Reid 标准 型 是 这 样 定义 的 ， 它 所 包含 的 每 个 方程 可 表示 成 如 下 的 形式 


3 Au=0 fy C— D Xf. 


D,V? = f; (R) 


21 


RE: 分数 微分 方程 和 非 线性 发 展 方程 中 问题 的 研究 


其 中 等 式 左边 为 首 导数 ， 右 边 为 关于 > RV 的 导数 的 函数 ， 且 满足 如 下 条 件 ， 
(1) 在 每 一 个 方程 R) 中 ， 右 边 JT 所 包含 的 导数 均 低 于 其 首 导数 DV; 
(2) 对 于 不 同方 程 中 出 现 的 首 导数 是 不 同 的 ; 
(3) 首 导数 的 任何 导数 均 不 出 现在 其 它 任何 方程 的 函数 中 ， 
利用 该 标准 型 的 思想 ， Reid 成 功 地 将 线性 PDEs 约 化 为 标准 型 931. 
定理 321 设 X 为 n 维 线性 空间 ， 4,C,D 为 XX 上 的 线性 算 子 , 若 C 
D 对 Vidal: 证 明 :， 定义 一 个 序 ; 设 X = {(21,20,...,2n)|ai € Ri = 1..,n},0 = 
(v, ..., v7), = (ul, u?, ...,u*) 
(1) 对 不 同 的 变量 u,v, 按 v > u 排序 ; 
(2) 对 同一 个 变量 u 或 v, 不 妨 用 v 来 定义 ， 若 满足 下 列 条 件 : 
(i) ord(a) > ord(D); 
(ii) ord(a) = ord(8) H i > j; 
(iii) ord(a) = ord(B),i = j, H a -6 的 第 一 个 非 零 分 量 为 正 ， 
则 Dart > Dgv, 其 中 
porn 


Or] Org? Orat 
(A) Æ Au = 0 无 解 ， 即 不 能 用 Taylor PARRA, RAE Au = 0 为 


一 矛盾 方程 组 ， 又 已 知 hu = 0 为 e 的 形式 相 容 条 件 的 一 部 分 ， 即 Au = 0 为 


a = (Qi, On) € N”, B = (Bis B.) € NF, Dav? = ord(a) = al 十 … 十 an > 0. 


v=0 
a 的 划一 等 从 方程 组 的 一 部 分 ， 因 此 该 等 人 方程 组 无 解 ， 进 而 | ME: 
f&. 所 以 CKerD= 9 5 KerA=®. 
(B) Xi Au =0 有 解 
在 前 面 所 定义 的 序 下 ， 因 为 C.D 为 线性 算 子 ， 根 据 Reid 的 标准 型 算法 ， 我 们 可 将 
Cv=u 
ih ; 约 化 为 一 个 等 价 的 标准 型 方程 组 
v= 
= filu, v), i-21,2,..,8, 
Davi = gilu), j=1,2,..,l. 


Ep 
G =u Dv = filu v), 1= 1,2,...,5, 
= H 
Dy v —gjí(u), j21,2,..L. 
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(B1) i Au =0 A5 Dyw = glu) 一 样 ， 则 


Cu 一 化 Dat fuv), i= 1,2,...5, 
e 
Dv 20 Au = 0, 


所 以 Au = 0 的 解 可 以 用 Taylor FERRARA, TURTLES v 也 存在 形式 将 级 数 解 ， 由 


上 式 知 ， 对 ve Ker4, au 使 得 | Cv=4 即 CKerDD KerA. (82) E hu=0 形 式 上 与 
D 


v=0 
Dyvi = glu) 不 一 样 , 因为 A 为 线性 算 子 , 则 根据 前 面 的 规定 的 序 , 总 能 将 Au = 0 约 化 为 
Cv -u 


一 个 标准 型 ， Deut = hu (u), Jl] Au —0 e Daut = h,(u) BR Au = 0 为 人 的 形 
v=0 


Cv -u 


式 相 容 性 条 件 , 因此 Daut = hilu) 也 是 | 的 形式 相 容 性 条 件 中 仅 含 4 的 部 分 . 


Dv = 
Das! = filu,v), t= 1,2,...,8, 


Cv-u 
; BR Dy) = g;(u) 
: Dyv = gilu), j—1,2,..l. 


又 已 知 在 同一 序 下 | 等 价 于 | 


v= 


为 oa 的 形式 相 容 性 条 件 中 仅 含 u 的 可 积 条 件 ， 根 据 Reid 的 算法 思想 ， 对 同一 


方程 组 在 同一 序 下 得 到 的 相 容 性 条 件 应 是 唯一 的 . 因此 Deut = helu) 会 Divi = glu) , 
进而 有 Au =0 € Dv = 9;(w). 这 又 归结 为 情况 (B1), 所 以 有 即 CKerD > KerA. 定理 得 
证 . 对 于 任意 给 定 的 一 个 微分 代数 方程 组 DPS = 0, 如 果 不 是 显然 矛盾 的 ， 从 偏 微分 方程 
的 理论 来 说 ， 很 难 判定 它 是 否 有 解 ， 从 微分 代数 的 角度 出 发 ， 可 以 通过 Rosenfeld 引 理 ， 
判定 它 所 对 应 的 代数 方程 组 是 否 有 解 ， 而 判定 代数 方程 组 是 否 有 解 ， 我 们 可 以 用 Grobner 
基 的 方法 ， 一 个 代数 方程 组 无 解 当 且 仅 当 它 的 Grobner 基 为 1 . 
R= Kuta e untu uh E m 个 微分 算 子 的 微分 环 , 4u, Cv, Du C 

R, 

HE: Au 是 指 : 在 它 的 每 一 个 微分 多 项 式 中 ， 只 出 现 uu ,ut 和 它 的 有 限 阶 导 
A. 

Dv 是 指 ， 在 它 的 每 一 个 微分 多 项 式 中 ， 只 出 现 vuv ovs 和 它 的 有 限 阶 导数 
u = Cu 是 指 一 个 变换 wi = Oiv,i = 1,2,.… t 

设 hu - 0 是 我 们 所 要 求解 的 方程 组 ,经 过 变换 u= Cv ， 得 到 目标 方程 组 Dv = 0 . 
WMR Dv = 0 有 和 解 ， 一 般 地 有 Cd 一 zero(Dv) C d- zero(Au) .车 无 解 ， 则 说 明 Au = 0 没 
有 w= Cv 类 型 的 解 ， 下 面 的 定理 给 出 了 二 者 相等 的 充 要 条 件 ， 

定理 3.2.2: 设 Au =0,u = Cu, Dv =0 , MAB Du, Co- AEA u 的 可 积 条 件 恰 为 
Au 的 特征 集 当 且 仅 当 d - zero(Au) = Cd- zero(Dv) , 这 里 恰 为 Au 的 特征 集 是 指 ，4vw 的 
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特征 集 与 系统 Dv, Cv-u 的 特征 集中 只 含 4 的 部 分 , 或 者 相同 , 或 者 可 以 互相 约 化 为 0， 
这 里 的 序 为 ， 在 前 面 规定 的 序 的 基础 上 ， 再 规定 ww < u< weno <w. 

WERS: 由 Rosenfeld-Grobner 算法 ， 首 先 判 定 如 =0，w= Cv，Dv =0 是 否 有 解 ， 
车 无 解 ， 则 说 明 4u = 0 没有 = Co 类 型 的 解 .下 面 假设 其 有 解 。 

(= ) 反 证 . i Au 的 特征 集 为 DCS. 假设 Dv, Cu-u 还 有 其 它 只 含有 u 
的 可 积 条 件 ， 记 为 Pu, ERER DCS 约 化 为 0, 设 Pu = Pu mod[DCS|. ,此 时 ， 系 
H Dv, Cv -vv 的 特征 集中 只 含有 u 的 可 积 条 件 正 好 是 DCS, Piu 的 特征 集 ， 记 为 
DCS,. ABRAM, d-—zero(DCS1) C d - zero(DCS) ， 且 是 真 包含 ， 否 则 Piu 就 会 被 
DCS 约 化 为 0 ， 这 说 明 使 得 Dv = 0 Cv =u 有 和 解 的 4 的 集合 为 4 — zero(DCS1) ,而 
d- zero(DCS) =d—zero(Au) , FA. 

(=) 由 Rosenfeld-Grobner 算法 , 可知 Au=0, u=Cv, Dv=0 AR. 设 它 的 特征 
Sy DCS, 而 d— zero(Au,u - Cv, Dv) = d- zero(DCS) , BERRIRA Dv, Cv-u 
的 只 含有 u 的 可 积 条 件 为 恰 为 Au 的 特征 集 ， 可 知 ， 使 得 Dv =0，Cv 一 u =0 有 和 解 的 4 
AREER HA d- zero(Au) . Bl Vu € d — zero(Au), 3v € d— zero( Dv), BEE u-Cv, HI [al 
理 可 证 反 包 含 关 系 成 立 ， 定 理 得 证 ， 

对 于 线性 变 系数 的 情形 ， 如 果 我 们 考虑 其 正则 解 的 情形 ， 也 就 是 说 ， 限 制 在 某 个 代 
数 流 形 上 ， 其 结果 和 上 面 的 结果 类 似 . 

对 于 非 线性 的 情形 , 如 果 我 们 在 局 部 意义 下 考虑 其 正则 解 的 情况 , 我 们 有 以 下 结论 . 

定理 3.2.3: 设 Au = 0.4 = Cv, Dv = 0, 如 果 系 统 Du ， Cu -% 的 只 含有 的 可 积 条 
件 能 被 Au 的 特征 集约 化 为 0， 则 d- zero(Au/J1) € Cd — zero(Dv/J2) ,其 中 Ji. Ja 231] 
为 hu 的 特征 集 和 Dv, Co — u 的 特征 集 (RE v 的 部 分 ) 的 15- FER (证 明 从 赂 ) . 

一 般 情 形 下 ， 在 非 线 性 偏 微分 方程 精确 求解 的 过 程 中 ， 构 造 其 C- D MARAE 
有 两 种 : 一 是 直接 法 ， 即 直接 从 方程 本 身 出 发 ， 直 接 构造 满足 需要 的 C- DX] 二 是 假设 
法 , 即 基于 一 定 的 需要 , 对 CD 对 的 构造 提出 一 定 的 假设 , 在 假设 的 基础 上 构造 需要 的 
C-D} 这 种 方法 目前 很 多 ， 例 如 精确 求解 中 的 . Lax Xt, Darboux BH, Backlund Æ 
换 、 Painleve 检 验 、 齐 次 平衡 法 、 相 似 约 化 等 等 . 运用 这 一 思想 , 解决 了 一 类 非 线 性 系统 
的 求解 问题 , 使 其 在 弹性 力学 、 电 动力 学 、 流 体力 学 、 量 子 力学 和 孤立 子 理论 中 得 到 了 广 
ZAH, 尤其 是 在 非 线性 系统 的 求解 中 , 运用 数学 机 械 化 的 方法 , 使 得 部 分 C — D XIRS 
构造 得 以 机 械 化 的 实现 , 大 大 节约 了 人 的 脑力 劳动 , 丰富 和 发 展 了 “ AC = BD” HRR. 


83.3 ”一 种 构造 的 Darboux 变换 及 其 应 用 


Backlund 变换 和 Darboux 变换 
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1885 年 ， 瑞 典 几何 学 家 Backlund 在 研究 复 常数 曲面 时 ， 得 到 了 Sine-Gordon 方程 的 


有 趣 的 性 质 . 
d u Alu! 是 Sine-Gordon 方程 wy = sinu 的 两 个 解 ， 它 们 之 间 有 如 下 的 关系 : 


sin(— >u 
2 


g 


这 就 是 Sine-Gordon 方程 的 Backlund 变换 .该 变换 给 出 了 从 Sine-Gordon 方程 的 一 个 解 得 
到 另 一 个 解 的 方法 9 多 ， 另 外 得 到 了 一 个 非 线 性 琶 加 公式 ， 


fy F Hh, an( uz Jl + ug. (3.3.2) 


ul 
hun = -un + ). (3.3.1) 


NN: u+ 
ug = Ug 28 sin( j 


ui 一 incen a. 


在 今天 看 来 , 这 个 公式 在 非 线性 理论 中 具有 重要 作用 , 但 在 当时 ,由 于 没有 得 到 广泛 的 应 
用 ,而 被 冷落 了 近 百 年 . 在 直到 20 世纪 60 年代， 由 于 非 线 性 光学 ， 晶体 位 错 等 许多 领域 
的 研究 都 与 Sine-Gordon 方程 有 关 . 这 时 Backlund 变换 重新 引起 了 人 们 的 注意 ,并 成 为 制 
作 多 个 孤立 子 的 重要 手段 ， 尤 其 令 人 瞩目 的 是 ， Backlund 变换 的 存在 ， 可 换 性 定理 各 非 
线性 释 加 公式 等 事项 不 是 Sine-Gordon 方程 所 专 有 的 。 1973 年 ， Wahiquist 和 Estabrook 
RR Kav 方程 P9) 

Uu, + utu + Ure = 0. (3.3.3) 


也 具有 Backlund $H, $ w= i ° udz, 它 满足 
Wi + 3u2 + Weer = 0. | (3.3.4) 


作 方 程 组 
wh = f -w,-i(w-wy, (3.3.5) 


wj = -w + (W — w')(Wee — wig) - 2(w2 + we + w? 

容易 验证 ， 当 u 满足 KdV 方程 时 ， 方 程 335) 完全 可 积 ， 其 解 w' 的 导数 = wv) 也 是 
(3.3.3) 的 解 . 不 仅 如 此 ， 对 KdV 方程 也 有 类 似 的 可 换 性 定理 ， 其 非 线 性 秋 加 公式 为 : 

wiz = wot iA - a (3.3.6) 
1976 Æ Wahlquist 和 Estabook!5] 提出 了 非 线性 方程 的 Backlund 变换 的 延 拓 性 ， 把 Back- 
lund 变换 , FER, 反 散 射 变换 统一 在 一 个 新 概念 中 - 拟 位 势 中 ，1983 年 ，Weiss Tabor 
和 Carnval 推广 了 常 微分 方程 的 Painleve 可 积 的 判定 方法 ， 提 出 了 偏 微分 方程 的 Painleve 
可 积 的 判定 方法 ， 并 用 其 来 获得 可 积 方程 的 Backund BR, UR Darboux FR. 在 一 百 
多 年 前 达 布 (J. G. Darboux) RAR, BEBE 


Pz + ud = rAd, (3.3.7) 
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在 达 布 变换 
$ = br + 0¢,% = u— 2, = u + 2[Inf (zÀ3i)]as, (3.3.8) 
og = fa 20 for + uf =f (3.3.9) 
no 
下 是 不 变 的 . 即 $ 满足 与 (3.3.7) 形式 相同 的 方程 
Ger + (u)ó = A(G). (3.3.10) 


Darboux 变换 的 基本 思想 为 ， 利 用 非 线性 方程 的 一 个 解 及 其 Lex 对 的 解 ， 用 代数 算法 及 
微分 运算 来 获得 非 线性 方程 的 新 解 和 Lax 对 相应 的 解 . 有 时 人 们 将 Darboux 变换 也 称 
为 Backlund $K, RAMAK Backlund 变换 的 Darboux 方法 . 关于 Backlund 变换 的 早 
期 工作 可 参考 文献 96. 97,98) 1975 Æ, Wadati 等 人 将 Darboux 变换 推广 到 mKdV 和 
Sine-Gordon 方程 991. 1986 年 ， 中 科 院 院士 谷 超 豪 等 人 将 Darboux 变换 推广 到 KdV 族 ， 
ANKS 族 及 (1+2)- 维 ， 高 维 方程 组 ， 并 且 将 Darboux 变换 应 用 到 微分 几何 中 的 曲面 论 和 
调和 映照 中 ， 另 外 ， 延 拓 法 及 局 部 高 阶 切 从 法 等 也 能 获得 Backlund 变换 (98, 100, 101] F 
换 性 定理 和 非 线性 释 加 公式 则 给 出 了 解 之 间 的 代数 运算 ， 可 由 已 知 解 求 得 新 解 。 但 实际 
上 越 往 下 计算 越 复杂 ， 这 就 限制 了 Backlund 变换 的 应 用 ， 而 不 变 定 理 和 非 线 性 秋 加 公式 
(89, 95) 则 给 了 解 之 间 的 代数 运算 可 由 已 知 解 得 到 新 解 。 胡 星 标 教授 在 这 方面 作 了 深入 的 
研究 工作 2021-1108) | AAT, Backlund 变换 已 成 为 研究 非 线 性 方程 的 有 力 工具 . 特别 , 3X 
于 有 限 维 可 积 系统 的 Backlund 变换 又 引起 了 人 们 的 广泛 重视 ,我 国 数学 家 曾 云 波 教 授 作 
了 许多 的 研究 工作 107 . 

王 明 亮 教授 和 李 志 斌 教授 提出 了 求 Backlund 变换 的 有 效 而 简便 的 方法 . 范 恩 贵 教授 
ARIAL, AREI, RRA. SVB SET FS LH (SLRABR 
论文 20) . 

最 近 ， 楼 和 胡 008 通过 对 Kadomtsev-Petviashvit 方程 进行 对 称 约束 得 到 了 0-1 
维 ，(2+1)- 维和 (3+1)- 维 的 可 积 的 孤立 子 方程 ,其 中 (2+1)- 维 高 阶 Broer-Kaup(HBK) 系 


统 是 ， 
Hg + (Has + H? - HH, + 3HW),, + 12(HG) 22 = 0, 


G;--A(Gz, --3H?G + 3HG, + 3GW)s —0, — (3.3.11) 
Gz - W, — 0, 
(1+1)- 维 高 阶 Broer-Kaup(HBK) 系统 是 : 
H, 4(Hrs + H? — 3HH, 6HG), = 0, 
Gi +4(Grr + 3H?G + 3HG, + 3G?), = 0. 


r= y, A (3.3.11) 可 以 约 化 为 (3.3.12) . 我们 知道 许多 学 者 已 经 研究 了 (2+1)- 维 
高 阶 Broer-Kaup(HBK) 系统 的 Peinlevé 性 质 ， 一 致 局 部 结构 和 精确 解 等 各 方面 ， 并 且 ， 
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(1-+1)- 维 HBK 系统 也 吸引 了 越 来 越 多 数学 工作 者 的 注意 力 ， 它 在 其 他 学 科 有 很 多 重要 
的 应 用 ， 近 来 ， 范 和 黄 等 (109, 110) 利用 不 同 的 Darboux 变换 得 到 (1--1)- 维 高 阶 Broer- 
Kaup(HBK) 系统 的 精确 解 . 因为 Darboux 变换 是 求解 偏 微 分 方程 的 最 有 效 的 方法 之 一 . 
因此 我 们 也 试图 利用 相应 谱 问题 规范 变换 ， 构 造 一 些 新 的 Darboux 变换 来 求解 (1+1)- 维 
高 阶 Broer-Kaup(HBK) 系统 ， 通过 严格 的 代数 算法 获得 具有 一 定 物理 意义 的 显 性 解 , 
上 述 (1+1)- 维 高 阶 Broer-Kaup(HBK) 系统 的 Lax 对 可 以 由 其 相应 的 谱 问 题 (是 谱 


参数 ) : 
1- H) -G 
3 2 T a 2 
oz = Mo, $ = ($1,032) , M ( 1 VM 
和 辅助 谱 问 题 
Nu Ny 
=N 1 一 s 
" g~ (m bd 

来 构造 ， 其 中 


Ny = -和 3 - GA 十 就 Haoz - 2Gz + H? + 2HG - 3H H3), 


Ni = GA? + (Ga + HG)À + ((Gs + HG); + H(Ge + HG) + 2G?], 


Na = -\? ~ HÀ + (Hz - 2G — H?), 


很 显然 ， Lax Xt (3.3.13),(3.3.14) 可 以 通过 作 规 范 变换 


转化 为 : 
Gs = Mó, M=(T,+TM)T™, 


和 
h= Nó, NS (Dn 4 TN)T-!, 


由 相 容 条 件 drt = dis 可 以 产生 : 


M, — Ns + [M, Ñ] = T(M; - N,+([M,N])T™, 


(3.3.13) 


(3.3.14) 


(3.3.15) 


(3.3.16) 


(3.3.17) 


(3.3.18) 


它 可 以 告诉 我 们 要 使 HBK 系统 (3.3.13) 在 规范 变换 下 解 不 变 , 最 重要 的 就 是 寻找 与 M,N 
有 相同 形式 的 矩阵 MÑ. HE, M,N "PÉSIBISSE HO G 也 映射 成 了 MÄ 中 的 五 


ma. 
我 们 构造 了 如 下 Darboux 变换 : 


a(Ata) œ 


dai be — (A +4) 
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EH a = (c 1)/2 和 5= (ec+D-Vz(c 天 -了 ME. a,b,c,d 是 由 如 下 的 代数 系统 决定 : 


alà + a)dy + abd = 0, 
alà + a) + abya = 0, 
beg; + (Ai  d)és = 0, 
des + ó(Ao + d)» = 0, 


(3.3.20) 


其 中 $= (1, $2)? = (da, ba)? 是 分 别 关 于 谱 参 数 和 ,和 2, 和 1 A An KREE (3.3.13), (3.3.13) 


的 两 个 基本 解 ， 
通过 直接 计算 我 们 容易 得 到 : 
PI A262V1 一 A191/o 
EET EE 


b 
e Q2 -Agaa 
E E 


_ (M1 — Ahy 

= 

d= Adaya — Adie 
A , 


A = diio ~ payi. 
HAA (3.3.19) 说 明 T(A) 是 关于 参数 和 的 二 次 多 项 式 ， 且 


detT(A) = (A +a)(A + d) — be. 


把 wbed 代 入 dcir(A) ,我 们 可 以 得 到 ， detT( 和 4) = 0 = 1,2). 


注意 到 detT (A) 的 首 系 数 为 1, detT (4) 又 可 以 改写 为 : 
detT(A) = (A — Aj)(À ~ A). 


从 上 述 事 实 我 们 可 以 得 到 以 下 的 命题 : 


命题 1. 在 (3.3.16) "HB AERE M = (T; +TM)T 和 和 矩阵 有 相同 的 形式 ， 即 : 


a -G 
1 -jo-H) 


其 中 
A-H- 


G = (e + 1)(b +G). 


(3.3.21) 


(3.3.22) 
(3.3.23) 


(3.3.24) 
(3.3.25) 


(3.3.26) 


(3.3.27) 


(3.3.28) 


(3.3.29) 


变换 (3.3.15) 和 (3.3.19):(9, H, G) — (9, H,G) 就 是 谱 问 题 (3.3.13) 的 Dauboux 变换 ， 


证 明 : 令 T = T*/detT, ffl 


(T5; - TM)T* = ( f11(X) fia) ) 


falà) falà) 
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其 中 7* Ree T 的 伴随 矩阵 ,容易 知道 mA) 和 fa) 是 关于 参数 和 的 3 次 多 项 
X fi) 和 fa 0) 是 关于 参数 和 的 2 次 多 项 式 . 我 们 容易 证 明 A A 是 flA j = 1,2) 
根 ， 因 此 deT(A) | fu09(53 = 52). 

基于 以 上 的 事实 ， 我 们 可 以 假设 : 


(T, + TM)T* = detT P(A), (3.3.30) 


其 中 a) (0) (0) 
P PY à t Py p12 ) 
o) ( pl pat py) 
Bopp ,3 =1,2k =0,1) 是 与 和 无 关 的 待定 函数 ， 则 方程 (3.3.30) 可 以 改写 成 


(T, + TM)T = PQJT. (3.3.31) 


比较 方程 两 边 入 ,项 的 系数 ， 可 以 得 到 : 


Pu = 3,P22 二 一 二 

pi 一 一 apii 十 二 一 各 一 等 = 一 后 
iD = Pol) tee pa 
pO -$p aG _ a me 
py =~ Spy te tha 


JAK (3.3.16) , (3.31) 可 以 看 出 : M= PA). 
命题 2. 在 Darboux BH (3.3.15), (3.3.19) F, WM a,b, ed 满足 方程 组 (3.3.21)-(3.3.25) , 
则 


G 
b, = Ga - —d- Hb (3.3.32) 
C; = (c 1)a + Hc- d (3.3.33) 
e = 8(c + D — h) + (G — G)a + T ;(6« + HG) + Hi] (3.3.34) 
成 立 ， 其 中 : 
= E 3 e 
le Uta - 16s E848 - SI aan 
h = 4(Hzs — 2G, H5 + 2HG — 3H H3). 
证 明 ， 从 关系 式 ， 
dis = $03 一 互 ) 轴 -Gdb2， 
edi e 
Qar 二 内 3(A1 H)óa, (3.3.36) 


Wiz = 302 — H)ji - GV», 
Vas = 91 — $02 — H)vs. 
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容易 计算 出 : 
(bi)z = $00 十 Dan- ~ Hory — G(gav + $12), 
(22)z = —1Qu + A2)ó2Ua + Hery + ($21 + $12), 
(页 如) = b(A — M2) pia — Gorp + div, (3.3.37) 
(dats )e = $(7À + Az) oot — Gap + Git, 
= (diva — 21); = 5 (Ar — As) (Vna + 42a). 


利用 方程 组 (3.3.21)-(3.3.25) ， 我 们 可 以 得 到 : 


s cur + 和) 和 而 — Host — Glive + dav] 


az = Andry 508 — Ao (se doi) 


EIS 39) ade + $26] - HO Gig] - M9 ossia - Jodi) 
- pat Aio sAn Meter _ ge 9 SM 


E^ Oa = Aa oad dada) 


=Ga~ © a - Hc 
c+l 
(3.3.38) 


用 同样 的 方法 和 技巧 ， 我 们 可 以 证 明 ct = (c+l)a+He-d, 此 外 ， 我 们 可 以 计算 出 在 辅 
助 的 谱 问题 中 的 【ga ji (0292): (的 oa) (0201), . SAT, ATER N 的 复杂 性 ， 这 个 过 
程 将 十 分 复杂 , 利用 计算 机 符号 计算 软件 Maple ， 我 们 可 以 证 明 (3.3.34) 从 而 完成 命题 证 
88. 

命题 3. XE DT 变换 (3.3.15),(3.3.19) F, (3.3.17) 中 的 矩阵 N SE N 有 相同 的 形式 ， 


Rp: T 
N=4 ( Nu Na (3.3.39) 


其 中 
Ny = -443 — GA + 4(Aiz ~ 2G, + Ë? + 2HG - 3H.) 
My = Gd? + (Gz + AG)A+ (Gr + HG), + B(G, + HG) + 264], 
Ny = -A* — Hd + (H, — 2G — H*), 
同时 原 函 数 豆 和 GHAR AOA G( 用 命题 1 中 的 所 用 的 方法 , 我 们 可 以 证 明 该 命题 ) . 
Darboux 变换 的 应 用 
接 下 来 ， 我 们 将 应 用 DT 变换 构造 HBK 系统 (3.3.12) 的 显 性 解 ， 它 们 具有 一 定 的 物 
_ 理 意义 . 为 了 简单 起 见 ， 我 们 令 万 = 0,G = 1 ， 选 择 三 个 基本 解 ， 


à; = cosh[£i(z,t)], $2 = ay sinhfEr(z, t)] + $Acosnlé (z, t)], 
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vy = sinh[£2(r,t), Vo = agcosh[£(z, t)] + 3Aasinh[£o(z. t), 

$3 = cosh|é3(z,t)], ba assinh(£a(z, t)) + pAscosh[Es(z, t)]. 
其 中 = -HA - 1)2, &(z,t) = -ailt + &t)(B = —4(02 + 2), (i = 1,2,3)), BI Ë = H — $ 
fi G=(c+1)(b+G). 
分 析 可 知 ， 如 果 Ay > Ar > Labod RA 由 式 (3.3.21)-(3.3.25) 决定 ， 则 (A(z, t), G(s t)] 
是 两 个 追赶 碰撞 孤立 波 ; 而 且 ， 如 果 A > 1,r2 < =i, a, b, cd , A = $104 $203 . 得 到 的 
新 的 解 [A.C] 是 两 个 迎面 碰撞 孤立 波 . 
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余 寄 ;分 数 微分 方程 和 非 线性 发 展 方程 中 问题 的 研究 
结 i 


力学 、 物理、 化 学 、 生 物 、 金 融 学 领域 所 引出 的 分 数 微分 方程 和 非 线性 发 展 方程 具有 
明确 的 实际 意义 ， 而 求解 这 些 方程 往往 是 比较 困难 的 .试图 解决 其 中 的 一 些 问 题 ， 本 文 
得 到 了 如 下 的 结果 和 一 些 仍 有 待 解决 的 想法 : 

一 、 构 造 了 一 种 新 的 Mitteg-Lefer 函数 ， 并 且 研 究 了 它 的 量 要 的 积分 变换 ， 包 括 
Laplace 变换 ， Mellin 和 逆 Fourier 变换 ， 为 求解 带 有 分 数 时 间 和 空间 导数 的 分 数 扩散 
- 波动 方程 握 供 了 较 有 效 的 工具 ， 

二 、 利 用 分 数 微分 形式 推广 了 向 量 分 析 中 的 基本 工具 一 - 梯度 ， 旋 度 和 散 度 到 所 谓 
的 分 数 梯度 ， 旋 度 和 散 度 . 

三 、 以 机 械 化 的 方法 为 指导 构造 了 一 种 Darboux 变换 ， 将 一 类 微分 方程 转化 为 代数 
方程 来 求解 ， 并 得 到 了 一 些 有 意义 的 新 解 . 

近来 出 现 了 越 来 越 多 的 分 数 发 展 方程 模型 (主要 是 线性 的 ) , 我 们 考虑 是 否 可 以 将 比较 
成 熟 的 非 线性 发 展 方程 的 理论 应 用 于 发 展 不 够 成 熟 的 分 数 微 分 方程 的 理论 ? 分 数 Maxwell 
方程 以 及 一 些 相关 方程 的 出 现 让 我 们 思考 ， 既 然 外 微分 形式 可 以 改写 电磁 学 的 经 典 方程 
组 ， 分 数 微分 形式 是 否 可 以 改写 这 些 新 出 现 的 具有 物理 意义 的 分 数 微分 方程 ? 这 些 分 数 
微分 算 子 在 什么 样 的 条 件 下 可 以 代数 化 , 机 械 化 ?.. .， 因此， 我们 的 研究 还 有 较 大 的 空 
间 . 
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